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In support of the Barents Sea Polar Front Experiment (BSPFEX) 
in September 1992 (Barents Sea Polar Front Group, 1992), the 
planned 224 Hz tomography signal transmissions from a near bottom 
sound source to a vertical receiver array consisting of 16 
hydrophones were simulated. Acoustic rays were traced to the 
receiver array at a range of 50 km using the NOAA Hamiltonian 
Raytracing Program for the Ocean (HARPO). Input to HARPO was a 
mathematical ocean environment based on historical bathymetric and 
sound speed data. Acoustic multipath arrival structure was 
constructed through eigenray searches and estimation of raytube 
spreading and surface and bottom losses. A resolvability analysis 
of the simulated arrival structure reveals that there are a total 
of 49 unique resolvable ray arrivals. Among them, 42 are from 
individual omnidirectional hydrophones and 7 from plane wave 
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ITI 


II. 


А. 


В 


Ե: 


III. 


о qA 


IV. 


A. 


B. 


C: 


© 


TABLE OF CONTENTS 


PHYSICAL OCEANOGRAPHY .... . 


INTRODUCTION е e e 9 9 9 9 9 9 
. WATER MASSES AND SURFACE WAVES 


BARENTS SEA POLAR FRONT . . . . 


ACOUSTIC EFFECTS 
INTRODUCTION 
ABSORPTION 
SPREADING AND REFRACTIVE EFFECTS 
SURFACE LOSS 


BOTTOM LOSS 


RAY THEORY ACOUSTICS (0.5. .Լ..92 
INTRODUCTION 
HAMILTONIAN RAY TRACING 


HARPO OVERVIEW 


а 


I. INTRODUCTION .. 7 1 ee 
A. OCEAN ACOUSTIC TOMOGRAPHY .. . 
В. BARENTS SEA POLAR FRONT EXPERIMENT 
C. THESIS OBJECTIVES AND APPROACHES 
D. THESIS OUTCINE EE Ն. ԽԵ 


° 


12 


15 
15 
15 
15 
16 


18 


20 
20 
20 


21 


D. MODELING THE BARENTS SEA ACOUSTIC ENVIRONMENT 
1. General 
2. Sound Speed Field, Bathymetry, and Sea 
Surface 


3 Summary “Go. . 


՞ 
՞ 
՞ 
՞ 


Մ. ARRIVAL STRUCTURE SIMULATION & ANALYSIS 
A. INTRODUCTION 
B. RAY PATH STRUCTURE . 
C ARRIVAL STRUCTURE 
D. ARRIVAL TIME STRUCTURE 
БВ, UALTIVE AMPLITUDE . . . . e è a œ 
F RESOLVABILITY ANALYSIS . . . 
Ն Ար Ե ԷԾՎԱՇԵԼՈ 
2. Method I - Individual Hydrophones 
a. Time Analysis 
D. o, Analysis 
c. Resolvable Rays 


3. Method II - Plane Wave Beamformer 
NI. CONCLUSIONS 
LIST OF REFERENCES 


APPENDIX A: LIST OF EXTERNAL PROGRAMS 


22 


22 


pe 


27 


29 
29 
29 
30 
32 
34 
39 
39 
40 
40 
41 
43 


46 


51 


53 


ЭЭ 


APPENDIX B: INDIVIDUAL HYDROPHONE RESOLVABILITY 
ANALYSIS TABLE 


APPENDIX C: LINE ARRAY RESOLVABILITY ANALYSIS TABLE 


INITIAL DISTRIBUTION LIST ОС СЕС 


Nr 


81 


92 


98 


LIST OF FIGURES 


Figure 1 Barents Sea bathymetry (after Eldhom and 
Talwami, 1977). Box centered at 74.6N 24.0E indicates 
the specific area of study. ΠΡ... 


Figure 2 Barents Sea Polar Front experiment configuration 
and the details of the bathymetry (20 m contour 
interval). ((A)-224 H2 Source, (B)-16 element vertical 
hydrophone array, (C,D)-400 Hz transceivers.) 


Figure 3 A three dimensional view of the study area. 
((A)-224 H2 source; (B)-Vertical ee ae - 
(C,D)-400 Hz transceivers.) . . шн, 


Figure 4 Geographical distribution of water masses in the 
Barents Sea (Loeng, 1991). The Barents Sea Polar 
Front is defined by the shaded line. Water masses are 
described in the text. Box centered at 74.6N 24.0E 
indicates the BSPFEX area. ο νο ο ee ο... 


Figure 5 Mean wave heights (ft) for the Barents Sea, Л 
- September (NAVOCEANO, 1990). : TES ἃ 


Figure 6 Temperature structure in °C (1°C contour 
interval) of the Polar Front in an area 100 km SW of 
the BSPFEX area. The shaded area is defined by the 3- 
4°C contour. (Johannessen and Foster, 1978). ὃ š 


Figure 7 Topographically controlled section of the 
Barents Sea Polar Front (broad line) around Bear 
Island (Johannessen and Foster, 1978) 


Figure 8 Surface loss at a sea state of 3 (Emblidge, 
1991) |  ^-—9 a 


Figure 9 Bottom loss along the 50 km path from the source 
location (A) to vertical hydrophone array location (B) 
(Kerr, (1990 and Blodgett, et al., 1987). : : 


Figure 10 Historical sound speed profiles. SSP1-North 


Atlantic Water, SSP2-Polar Бопе SSP3-Arctic 
Water.(NAVOCEANO MOODS, 1991). NON. WU VN 


vii 


Tor 


I2 


13 


14 


17 


1, 


24 


Figure 11 Linear interpolation of historical sound speed 
profiles. А1-А5 аге defined in the text. 
(Emblidge,1991). š »9պ.ծ..ծ....-.....--- 


Figure 12 Grid overlay for bathymetry model. 


Figure 13 Orientation of modeled sound speed profile. (P3 
- Arctic Water, P2 - Polar Front, Pl - North Atlantic 
Water, Al - A5 represent the linear interpolated 
averages.) © . ... ез. ο T o ЭО 


Figure 14 Modeled sound speed field along the track from 
Source А to hydrophone array B. m. αμ. . L. 


Figure 15 Ray path along a vertical slice from source (A) 
to receiver array (B) for a launch angle of 5.00?. 


Figure 16 Ray path along the track from source (A) to 
hydrophone array (B) for a launch angle of 8.89?. 


Figure 17 Ray path along the track from source (A) to 
hydrophone array (B) for a launch angle of 10.40°. 


Figure 18 Arrival depth versus launch angle for the near 
bottom sound source (A). M ле 


Figure 19 Arrival depth versus arrival angle for the near 
bottom sound source (A). . j. -...... 3 


Figure 20 Arrival depth versus arrival time for the near 
bottom sound source (A). 


Figure 21 Arrival angle versus arrival time for the near 
ԵԾԵՇԵԾՈ Տօսոմ Տօսօ6լ"')-ԶշռԴեզ-Շշռ՞շռշՒւպշպշՒշեե.եԼեԼքշԿ'խԽՒխԿ-օՄԽ/,Խ7 


Figure 22 Launch angle versus arrival time for the near 
bottom sound source (A). 


Figure 23 Individual ray transmission loss versus launch 
angle for a near bottom sound source (A) at a sea 
state of 3. © ез ЖОО Ое 


Figure 24 Individual ray transmission loss versus launch 
angle for a near bottom sound source (A) at a sea 
state of 5. ......աջ- 


Figure 25 Individual ray transmission loss versus arrival 


angle for a near bottom sound source (A) at a sea 
state of 3. ο առ աՆ... ԵԱ 


ү: 


24 


25 


26 


21 


30 


EM 


3d 


32 


33 


34 


35 


35 


37 


37 


38 


Figure 26 Relative amplitude of all rays from source (A) 
arriving at a vertical surface 50 km away at location 
(B) for a sea state of 3. 


Eure 27 SNR versus c, for two pulse widths. The 224 Hz 
source has a d of 62. 5 ms whereas the 400 Hz source 
has a d of 10 ms. 


Figure 28  Simulated arrival structure at ее #8 
(depth 220 m) at a sea state of 3. . . . . pdt 


Figure 29 Beam pattern for modeled linear array with beam 
Steered to 20 degrees off broadside. е с 


Figure 30 Simulated beamformed arrival structure in 


sector 13 (+13.5° grazing angle). R9 and R10 are 
resolvable and have o, < 6.3 ms. -————— 6 


12 


39 


43 


46 


47 


30 


LIST OF TABLES 


Table I LOCATIONS AND CHARACTERISTICS OF ACOUSTIC 
ELEMENTS . . . . 2M 


Table II BARENTS SEA WATER MASSES 


Table III INDIVIDUAL HYDROPHONE RESOLVABLE RAY 
ARRIVALS WITH o, < 6.3 ms 


Table IV PLANE WAVE BEAMFORMED RESOLVABLE ARRIVALS 
WITH Go, 6.3 ms .. eee oe 


10 


45 


49 


I. INTRODUCTION 


A. OCEAN ACOUSTIC TOMOGRAPHY 

Ocean tomography is an acoustical method to monitor the 
ocean interior (Munk and Wunsch, 1979). Like the computer 
assisted tomography (CAT) scans used in medicine and seismic 
tomography used in geology, ocean tomography employs beams of 
energy to create a three-dimensional image of the volume 
traversed. In CAT scans, these energy beams consist of X- 
rays; in seismic tomography, shock waves from earthquakes or 
explosions; and in ocean tomography, low frequency sound 
waves. As sound waves travel through the ocean they gather 
information on the changes in ocean temperature and currents. 
The data are in the form of changes in sound pulse travel 
times. Using inverse techniques the best estimate of the 
ocean structure is constructed. 

An ocean acoustic ους problem can be divided into 
two parts. The first is known as the "forward" and the second 
as the "inverse" (Munk and Wunsch, 1979). The forward problem 
establishes the physical relationship between data and the 
unknown ocean e cet Simulation studies using this 
relationship can be used to investigate signal and array 


design issues. BM NVersedproDlemwsdeals with the 


reconstruction of the unknown ocean structure based on the 
gathered data and the established forward relationship. 
The feasibility of using tomography for ocean monitoring 

depends on four factors (Munk and Wunsch, 1979): 

1. Stability 

2. Resolvability 

3. Identifiability 

4. Signal to Noise Ratio (SNR) 
Useful time series of acoustic travel time can only be derived 
from stable and resolvable multipath arrivals over successive 
transmissions. Identification (i.e. association) of the 
observed arrivals with modeled arrivals is required to 
determine the geometry of the acoustic paths. SNR determines 


the accuracy of the travel time measurement. 


B. BARENTS SEA POLAR FRONT EXPERIMENT 

The Barents Sea Polar Front Experiment (BSPFEX) is planned 
for August 1992 (Barents Sea Polar Front Group, 1992). The 
experiment will be conducted jointly by the Naval Postgraduate 
School (NPS), Woods Hole Oceanographic Institution (WHOI), and 
Science Applications International Corporation (SAIC). The 
experiment will take place within an approximately 70 x 80 km 
polygon centered at about 100 km east of Bear Island and 
Situated over the steep northwestern slope of Bear Island 


Trough as shown in Figure l. 


This field study will emphasize small-scale to mesoscale 
processes, and will utilize both traditional oceanographic 


measurements and acoustic tomographic techniques. 


20" 
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40° 45 





Figure 1 Barents Sea bathymetry (after Eldhom and Talwami, 
1977). Box centered at 74.6N 24.0E indicates the specific 
area of study. 


Figure 2 shows the experimental configuration and details of 
the bathymetry. The tomography system will consist of a 224 Hz 
source (A), a vertical array composed of 16 hydrophones (B), 
and two 400 Hz transceivers (C,D). Figure 3 shows a three 
dimensional view of the test area. Some of the 
characteristics of the sources and receivers are summarized in 
Table I. This study will focus on the 224 Hz source and the 
hydrophone array. The 224 Hz source is mounted 3 m off the 
bottom on an oceanographic mooring. The hydrophone array 
consists of 16 hydrophones mounted at a 10 m interelement 
spacing on a vertical oceanographic mooring with the top 
hydrophone at 150 m. 


Table I LOCATIONS AND CHARACTERISTICS OF ACOUSTIC ELEMENTS 





LAT LON FREQ BAND 
MOORING TYPE (deq N) (deq E) (ΒΖ) WIDTH(Hz) 
A Source 74.91 24.39 224 16 
B Hyd 
Array 74.51 2928 N/A N/A 
C Trans- 
| Ceiver 74.28 23:85 400 100 
р Тгапѕ- 
ceiver 74.67 22.95 400 100 











The scientific objectives of the tomography experiment can 


be summarized as follows (Barents Sea Polar Front Group, 


1992): 
1. Provide a detailed physical description of the front. 
2 Enhance the understanding of dynamics of the front, 


including frontogenesis and its influence on regional 





Figure 2 Barents Sea Polar Front experiment configuration and 
the details of the bathymetry (20 m contour interval). ((A)- 
224 Hz Source, (B)-16 element vertical hydrophone array, 
(C,D)-400 Hz transceivers.) 


oceanographic processes. 

3s Assess the ability of acoustic tomography to define 
frontal and associated mesoscale features. 

4. Provide improved dece m of shallow water 
acoustic propagation in this region and the effect of the 


environment on acoustic ASW operations. 


Depth (x100m) 





Figure 3 A three dimensional view of the study area. ((A)-224 
Hz source; (B)-Vertical hydrophone array; (C,D)-400 Hz 
transceivers.) 


C. THESIS OBJECTIVES AND APPROACHES 

The objective of this thesis was to address the 
tomographic forward problem in the Barents Sea. The study 
Simulated and examined the expected multipath arrival 
Structure of the 224 Hz signal at the hydrophone array. The 


Simulation used a three-dimensional ray theory acoustics 


model. The major focus was to investigate resolvability and 
ieS relation to SNR. 

In simulating the multipath arrival structure, the 
Semeuter program HARPO (Hamiltonian Acoustic  Raytracing 
Program for the Ocean) was used to trace a fan of rays to the 
hydrophone array location. Input to HARPO was a mathematical 
ocean based on a high resolution bathymetry chart (Norsk, 
1986) and historical sound speed profiles. Transmission loss 
for each ray was then computed using computer programs 
external to HARPO. Eigenrays and hence arrival structure were 
then determined for each hydrophone of the vertical array by 
analysis of the rays' arrival depths and times. 

Resolvability of the eigenray arrivals at the hydrophone 
array were examined for two cases. The first case considered 
each hydrophone of the vertical array as an independent 
omnidirectional receiver. The second case considered the 
hydrophones of the vertical array together as a plane wave 


beamformer. 


D. THESIS OUTLINE 

The remainder of this thesis consists of five chapters. 
Chapter II describes the physical oceanography of the Barents 
Sea including water masses, surface waves, and Polar Front 


features. Much of this chapter is taken from Emblidge (1991). 


Chapter III describes the acoustic properties of the 
Barents Seas. Here the method of calculating transmission 
loss of individual acoustic arrivals is discussed. 

In Chapter IV a review of ray theory, the ray tracing 
program HARPO, and the environmental models used аге 
presented. Chapter V provides the simulated arrival 
structure results, transmission loss calculations, апа 
resolvability analysis. 


Chapter VI discusses the conclusion of this study. 


II.  PHYSICAL OCEANOGRAPHY 


A. INTRODUCTION 

The Barents Sea is bordered to the south by the coasts of 
Scandinavia and Russian Republic, to the north by the Svalbard 
Archipelago and Franz Joseph Land in the southern edge of the 
Arctic Ocean. It is bounded on the eastern side by Novaya 
Zemlya. Its western boundary is nearly open and can be 
approximated by the 15 degree east meridian (Figure 1). With 
an average depth of 230 m and a maximum depth of 500 m, the 
Barents Sea is among the shallowest seas of the world ocean 
(Klenova, 1966). 

The  Barents Sea contains a complex  oceanographic 
structure. The meeting of Polar and Atlantic water masses to 
the east of Bear Island forms the Barents Sea Polar Front. 

The oceanographic conditions described in this chapter are 
those that are expected to exist at the time of the BSPFEX 


(August and September). 


B. WATER MASSES AND SURFACE WAVES 
Currents in the region transport water of both Arctic and 
Atlantic origin. These water masses have vastly different 


temperature and salinity characteristics (Table II). 


Table II BARENTS SEA WATER MASSES (after Loeng, 1991) 





TYPE T (deg C) S (psu) 
Arctic Water (AW) <0.0 34.2 - 34.8 
Polar Front Water (PW) -0.5 Со 2.0 34.8 - 35.0 


North Atlantic 
Water (NAW) >3.0 235.1 





Figure 4 shows Loeng's model of the geographic regions 
occupied by the major water masses. The abbreviations are AW 
for Arctic Water, SBW for Svalbard Bank Water, NAW for North 
Atlantic Water, BW for Bottom Water, BSW for Barents Sea 
Water, and CW for Coastal Water. The BSPFEX region is 
indicated by a small square in Figure 4. 

As shown in Table II the NAW flowing ο... introduces warm 
saline water while AW flowing south introduces cold, 
relatively fresh water into the Barents Sea. The confluence 
of these flows creates the Barents Sea Polar Front. As in 
many frontal situations in the ocean, the front has a complex 
horizontal and vertical structure. 

An important effect on sound transmission is the roughness 
of the sea surface. Many acoustic paths are expected to be 
surface interacting because of the shallowness of the sea. 
Figure 5 shows the mean wave height in feet, in the summer 
months of July - September (NAVOCEANO, 1990). Note that the 


BSPFEX area has mean summer wave heights of about 3 feet (1 





Figure 4 Geographical distribution of water masses in the 
Barents Sea (Loeng, 1991). The Barents Sea Polar Front is 
defined by the shaded line. Water masses are described in the 
text. Box centered at 74.6N 24.0E indicates the BSPFEX area. 


meter). The effects of sea surface roughness on the 


scattering of sound will be discussed in Chapter III. 
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Figure 5 Mean wave heights (ft) for the Barents Sea, July - 
September (NAVOCEANO, 1990). 


C. BARENTS SEA POLAR FRONT 

Across the western interior of the Barents Sea a front 
exists due to the adjacent position of the NAW to the south 
and AW to the north (Figure 4). The front near Bear Island 
has been studied in some detail by Johannessen and Foster 
(1978). Figure 6 shows the structure of the front, 
annotated by the shaded area, with a complex vertical and 
horizontal structure. 

Johannessen and Foster (1978) suggested that the mean 
position of the front was approximately locked to the outer 


part of the Svalbard Shelf by the 100-m isobath (Figure 7). 
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Figure 6 Temperature structure in °С (1?C contour interval) 
of the Polar Front in an area 100 km SW of the BSPFEX area. 
The shaded area is defined by the 3-4?C contour. (Johannessen 
and Foster, 1978). 


Figure 6 shows a baroclinic subsurface structure of the front 
and that the surface signature is nearly absent. 

The frontal core oscillates with the tides and seasons. 
The tidal oscillation is on the order of 10 km/cycle around 
Bear Island (Johannessen and Foster, 1978). The seasonal and 
yearly positions vary on the order of 50 km/year (NAVOCEANO, 
1991). 

The Barents Sea Polar Front is characterized by a change 
їп temperature of about 5°C and salinity of 1 psu over the 100 
km of its horizontal extent near Bear Island (Dickson et al., 


1970). 
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Figure 7 Topographically controlled section of the Barents 
Sea Polar Front (broad line) around Bear Island (Johannessen 
and Foster, 1978) 


The BSPFEX location is about 100 km NE of the area studied 
by Johannessen and Foster. Frontal conditions in the BSPFEX 
area are expected to be similar to those found by Johannessen 


and Foster. 
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III. ACOUSTIC EFFECTS 


A. INTRODUCTION 

Understanding the propagation of sound in the ocean 
requires a knowledge of the properties of the ocean medium, 
its boundaries and their influence on sound propagation. 
Among the many environmental factors the principal ones 
include: absorption, spreading and refraction, and 


surface/bottom interactions. 


B. ABSORPTION 

Absorption is the loss of acoustic energy along its path 
due to the conversion of acoustic energy into thermal energy. 
This absorption loss (TL,.) can be approximated by the 


following equation: 


TL pa = IY (1) 


where r is the ray path distance in meters and y is the 
absorption coefficient. For 224 Hz, y 185 approximately 
0.000006 dB/m (Kinsler, et al., 1982). So for a distance of 


50 km TE is 0.3 dB. 


S 


C. SPREADING AND REFRACTIVE EFFECTS 
Emesdsugosdndesbetractrve loss (ΠΤ) 1s the loss ог 


apparent gain of acoustic energy over the traveled path due to 


15 


the defocusing or focusing of adjacent rays as they track 
away from a small spherical source. TL,, was calculated using 
a FORTRAN computer algorithm 'cordat.f' (Appendix A) for each 
ray by comparing the separation between adjacently launched 
rays. | Assuming cylindrical symmetry, the equation for TL, 


from C.S. Chiu (1992) is 


E mg 
ТІ = 10109--хвусов бу (2) 


where r is the range between source and receiver, h is the ray 
tube cross sectional distance at the receiver (the cross 
section is orthogonal to the eigen-ray path), O0 is the launch 
angle of the ray and A0 is the difference in launch angle 


between adjacently launched rays. 


D. SURFACE LOSS 

Due to the high potential for acoustic ray interaction 
with the surface in a shallow ocean, the losses at the surface 
can be important. Acoustic energy is lost at each ray 
reflection with the surface due to the scatter of energy out 
'of the ray path. This loss in energy can be described by εν. 
reflection coefficient at the sea surface which is the ratio 
of ray energy reflected to incident ray energy. 

Using an equation in Clay and Medwin (1977), for 
Scattering loss by a rough surface having a gaussian 


distribution, Emblidge (1991) calculated the surface loss per 


ΤΕ 


bounce corresponding to a rms waveheight of 1 m or a sea state 


աթ Figure 8 is a plot of the calculation. 
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Figure 8 Surface loss at a sea state of 3 (Emblidge, 1991) 
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The total loss due to surface interactions, TL can be 


sfc' 
computed by summing the losses of each surface interaction, 


sloss,, that the ray encounters over its path: 


TL 


Sfc 


N 
= Y sloss; (3) 


151 


E. BOTTOM LOSS 

The acoustic energy lost by the ray interacting with the 
bottom is due to the transmission of some of the acoustic 
energy into the bottom and the scattering of acoustic energy 
out of the ray tube. The loss of energy can be described by 
the reflection coefficient at the bottom which is the ratio of 
acoustic energy reflected to that incident at the bottom. 

To get an estimate of the bottom loss, the U.S. Navy 
standard Bottom Loss Upgrade (BLUG) is used (Kerr, 1990). 
BLUG is an interactive routine to compute a bottom loss curve 
from user-defined geoacoustic inputs. Using the standard Navy 
BLUG computer algorithm and geoacoustic parameters from an 
arctic acoustic model called ICECAP (Blodgett, et. al., 1987) 
a bottom loss profile for the BSPFEX area was computed 
(Keenan, 1992). The results for 224 Hz are shown in Figure 9. 

As can be seen in Figure 9 bottom loss per bounce for 
grazing angles less than 25 degrees are less than 7 dB for 
ranges less than 18 km. The bottom loss per bounce is about 


0.01 dB for ranges greater than 18 km. 
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Figure 9 Bottom loss along the 50 km path from the source 
location (A) to vertical hydrophone array location (B) 
(Kerr, 1990 and Blodgett, et al., 1987). 


If a constant loss per bottom bounce is assumed then the 
iaesedue to the bottom for the entire ray path, TL,., can be 


expressed as: 


lives N(DIOSS) (4) 


where N is the number of bottom bounces the ray encounters and 


bloss is the transmission loss in Db per bounce. 
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IV. RAY THEORY ACOUSTICS 


A. INTRODUCTION 

The propagation of sound in an ocean medium can be 
represented as rays. Ray theory provides а visual 
representation of the paths taken by sound energy and 
graphically illustrates how various ocean parameters affect 
the path of each ray. The ocean acoustic model used in this 
: thesis is the Hamiltonian Acoustic Raytracing Program for the 
Ocean (HARPO). 

HARPO is a three-dimensional acoustic ray theory computer 
model. A significant advantage of HARPO is that it deals with 
continuous sound speed and bathymetric structure. This 
continuous treatment overcomes the problems of false caustics. 


(Jones et al., 1986) 


B. HAMILTONIAN RAY TRACING 

Acoustic energy travels in the form of compressional 
waves. In the high frequency approximation waves behave like 
particles. The path of these "acoustic particles" can then be 
determined by integrating a differential form of the wave 
. equation. Hamilton's equation, which governs changes of 
position and momentum in mechanical systems, are applicable 


to sound propagation at high frequency (Lighthill, 1978). 
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Hamilton's equation has the general form: 











ш i-1,2,3 

dt да; որ (5) 
а πο. 

dt Op, d 


where H(p,,Po,P3i9,,92/93) is a Hamiltonian function describing 
the total energy of a system in terms of generalized 
coordinate system p and momenta q. For acoustic application 
q is the wave number vector and p is a coordinate system. In 
HARPO the coordinate system is spherical polar. Solutions to 
equation (5) are obtained by choosing initial values for the 
six values of q, and p, and integrating this system of six 
differential equations. For sound propagation in the ocean 


the Hamiltonian takes on the form: 


прос) =0=с(р,) а (6) 


where o is the angular wave frequency, c(p,) is the sound 
speed field and q^ is the magnitude squared of the wave number 
vector (Jones, et al., 1986). Along the raypath, the 


Hamiltonian is defined to be zero. 


C. HARPO OVERVIEW 

HARPO is a computer ray tracing algorithm that numerically 
integrates Hamilton's equations. It requires, for input, a 
continuous three-dimensional representation of the sound speed 


field, and a continuous two-dimensional representation of the 
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upper and lower reflecting surfaces. The upper and lower 
surfaces are the ocean surface and bottom respectively. 
Gridded sound speed and bathymetry fields can be made 
continuous by the use of empirical orthogonal functions 
(EOF's) and splines using computer subroutines external to 
HARPO. These routines originated from Newhall et al. (1987). 
HARPO does not compute signal amplitude or eigenrays. 
Amplitude and eigenray determinations are made by external 
programs applied to the HARPO output files called DOUTP and 
RAYSET. The documentation of HARPO by Jones et al. (1986) 
provides a complete description of the mathematics and 


computer coding of HARPO. 


D. MODELING THE BARENTS SEA ACOUSTIC ENVIRONMENT 
l. General 

The accuracy with which HARPO calculates ray paths is 
primarily dependent on the accuracy with which the ocean is 
described by the input models. Chapter II described the 
environment of the BSPFEX area. Modeling the environment for 
the simulation of acoustic ray arrival structure was 
accomplished by selecting environmental data that provided 
adequate mathematical description for the sound speed field 


and bathymetry. 
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2. Sound Speed Field, Bathymetry, and Sea Surface 

The bathymetry input to HARPO is shown in Figure 12. 
The square box was manually gridded into 13 subdivisions along 
each side with a 7.5 km spacing. The bottom depths were then 
visually read off the chart (Norsk, 1986) at the gridded 
intersections. This gridded bathymetry input was splined to 
provide a continuous topography using FORTRAN computer 
subroutines "ЕДТ ՆՈ: (Emblidge, 1991), accgrEIdder.tf' 
(Newhall, et al., 1987), and 'bottom.f' (Newhall, et al., 
1987). Appendix A contains a copy of 'tgridin.f' which has 
been extensively modified for this research. Programs 
Begridder.f' ane 'bottom.f' only required a one line 
modification to change the size of the horizontal and vertical 
grid spacing. 

Figure 10 shows the three sound speed profiles (SSPs), 
extracted from the NAVOCEANO MOODS data base, chosen to 
represent ES on the North Atlantic side of the Polar 
Front (SSP1), front interior (SSP2), and on the Arctic side 
553), respectively. See Figure 4 for the water mass 
Orientation. From the three SSPs described above, five others 
were created by Emblidge (1991) by linear interpolation, as 
shown in Figure 11 (A1,A2,A3,A4,and A5), to provide a gradual 
change in sound speed on each side of the front. 

The three-dimensional sound speed field was generated by 
tying a specific SVP or linear interpolated SVP to a diagonal 


of the bathymetry field (Figure 12). Figure 13 shows the 
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Figure 10 Historical sound speed profiles.  SSPl-North 
Atlantic Water, SSP2-Polar Front, SSP3-Arctic Water.(NAVOCEANO 
MOODS, 1991). | 
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Figure 11 Linear interpolation of historical sound speed 
profiles. А1-А5 аге defined in the text. (Emblidge,1991). 
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Figure 12 Grid overlay for bathymetry model. 


orientation of the SVP's in the modeled BSPFEX area. This 
sound velocity input was then used to generate a continuous 
sound velocity field by EOF's and splines using FORTRAN 
Esuputer subroutines 'gridin.f' (Emblidge, 1991), 'gridder.f' 
(Newhall, et al., 1987), and 'cgrid.f' (Newhall, et al., 
1987). Appendix A contains a copy of 'gridin.f' which has 
been extensively modified for this work. Programs 'gridder.f' 
and 'cgrid.f' only required a one line modification to change 


the size of the horizontal and vertical grid spacing. This 
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Figure 13 Orientation of modeled sound speed profile. (P3 - 
Arctic Water, P2 - Polar Front, P1 - North Atlantic Water, Al 
- A5 represent the linear interpolated averages.) 


generated a modeled Polar Front about 40 km wide, oriented 
along the bathymetry contours, centered approximately between 
the source and receiver. A contour plot of the modeled SSP 
field along the track between the 224 Hz source (A) and the 
vertical array (B) is shown in Figure 14. This modeled cross 
front vertical slice compares favorably to actual data as 


shown in Figure 6. 
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Figure 14 Modeled sound speed field along the track from 
source A to hydrophone array B. 


ime ssea Surtace was modeled as a flat sphere concentric 
-with the earth with a radius of 6370 km. 

E Summary 

The environmental model designed provided for a 40 km 

wide front between the source (A) and the receiver array (B), 
with the source in the cold Arctic water, and the receiver 
array in the warm North Atlantic Water. Table I provides a 
description of the source/receiver positions. 

In the simulation, the source was placed about 3 m off the 


bottom in about 163 m-of water. The receiver was placed 50 km 


2 


away in about 318 m of water (Figure 3). Тһе track from the 


source to the receiver array has a bearing of 151.8? true. 
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V. ARRIVAL STRUCTURE SIMULATION & ANALYSIS 


A. INTRODUCTION 

To simulate the propagation of the 224 Hz pulse signal 
from a near-bottom sound source, through the front to the 
vertical receiver array, ray traces were performed for launch 
angles between 0.00 to 25.00° at 0.01° increments. The rays 
were terminated after passing 50 km in horizontal distance. 
Since HARPO does not stop the rays exactly at 50 km, a 
correction was made via the FORTRAN program ‘'cordat.f' to 
"back-up" the rays and compute arrival depth and time at the 


50 km receiver surface. 


B. RAY PATH STRUCTURE 

Rays with launch angles less than 8.89° are refracted 
bottom-reflected (RBR) as expected due to the downward 
refracting sound profile (Figure 10). Figure 15 shows a ray 
launched at 5.00°. The first ray to interact with the sea 
КО асос опе ταν launched at 8.89°. Figure 16 shows the ray 
launched at 8.89 degrees. The first ray to become completely 
surface-reflected and bottom-reflected (SRBR) is the ray 
launched at 10.40 degrees. Figure 17 shows the ray launched 


at this angle. 
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Figure 15 Ray path along a vertical slice from source (A) to 
receiver array (B) for a launch angle of 5.00?. 
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C. ^ ARRIVAL STRUCTURE 

For the 2500 rays traced between 0.00 and 25.00? various 
plots were made to study the arrival structure. The launch 
angle versus arrival depth plot is shown in Figure 18. The 
number of ray arrivals at any one receiver in the vertical 
array can be found by drawing a horizontal line across Figure 
18 at the corresponding receiver depth.  Eigenrays are given 


by the intersection of the horizontal line with the curve. 
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Figure 16 Ray path along the track from source (A) to 
hydrophone array (B) for a launch angle of 8.89?. 
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Figure 17 Ray path along the track from source (A) to 
hydrophone array (B) for a launch angle of 10.405. 
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Figure 18 Arrival depth versus launch angle for the near 
bottom sound source (A). 


Figure 18 shows more arrivals for receivers deeper than 150 m. 
The hydrophone array occupies the water column between 150 m 
and 300 m. 

The arrival angle versus arrival depth plot is shown in 
Figure 19. It is interesting to note that the rays tend to 
arrive at some preferred directions independent of arrival 


depth, especially for arrival angles larger than 10?. 


D. ARRIVAL TIME STRUCTURE 
Examination of the arrival times of the acoustic rays 
reveals information about temporal resolvability of acoustic 


signals for tomographic inversion. To accomplish this 
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Figure 19 Arrival depth versus arrival angle for the near 
bottom sound source (A). 


examination various plots were made of the 2500 rays launched 
between 0.00 and 25.00?. 

The arrival time versus arrival depth plot is shown in 
Figure 20 for all rays launched. This plot indicates that the 
later acoustic arrivals are associated with a more tilted wave 
front. The figure also indicates that the duration of the 
multipath arrival structure is about 4 s. 

The arrival time versus arrival angle plot is shown in 
ազատ 21. In this plot one can see that the earliest 


arrivals come in at small grazing angles. 
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Figure 20 Arrival depth versus arrival time for the near 
bottom sound source (A). 


The arrival time versus launch angle plot is shown in 
Figure 22. Consistent with Figure 21, this plot shows 


generally that rays with smallest launch angles arrive first. 


E. RELATIVE AMPLITUDE 

The ray amlitude of each ray that arrives at the vertical 
receiver surface 50 km away from the source can be computed 
using the transmission loss analysis discussed in Chapter III. 
The source level for the 224 Hz source (A) is 183 dB so the 
ray amplitude of the i'th ray, RAMP,, may be computed as 


follows: 
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Figure 21 Arrival angle versus arrival time for the near 
bottom sound source (A). 
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Figure 22 Launch angle versus arrival time for the near 
bottom sound source (A). 
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RAMP, = 183-TLeoe, (7) 


where TL, 1S the totalatransmicsion Wocceor πο μα του 
The total transmission loss for each ray is a sum of the 
energy losses due to absorption, spreading and refraction, and 
surface/bottom interactions as discussed in Chapter III. 
Neglecting absorption (which is insignificantly small) then 


the total transmission loss τος. ls 


ТІ. F IC | MM (8) 


Using Equations (2), (3), (4), and (8) the total 
transmission loss was calculated for each ray using 
'cordat.f'. For simplification the bottom loss was taken to 
be a constant of 0.5 dB per bounce (i.e., an average over the 
50 km track). Figures 23 and 24 are plots of the results of 
the transmission loss calculation for launch angles from 0 to 
25° at 0.01° steps over the 50 km horizontal distance for sea 
states 3 and 5 respectively. These figures show a strong 
dependence of acoustic energy loss due to interactions with 
the sea surface. Clearly, ray angle and sea state acoustic 
energy losses are quite significant for launch angles 213° at 
a sea state of 5. The corresponding mean wave heights is 4 m. 
Due to the large surface losses, the horizontal axis of Figure 
24 was intentionally terminated at 15?. 

Figure 25 is a plot of transmission loss versus arrival 


angle of all rays arriving at the vertical receiver surface 50 
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Figure 23 Individual ray transmission loss versus launch 
angle for a near bottom sound source (A) at a sea state of 
3. 
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Figure 24 Individual ray transmission loss versus launch 
angle for a near bottom sound source (A) at a sea state of 
5. 
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km away. It shows consistently that rays arriving at narrow 


angles will encounter the least transmission loss. 
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Figure 25 Individual ray transmission loss versus arrival 


angle for a near bottom sound source (A) at a sea state of 
Be: 


A plot of ray amplitude versus arrival time of all rays 
arriving at the vertical receiver surface at a sea state of 3 
is shown in Figure 26. This figure shows that the trend of 
the amplitude decreases with increasing ray arrival time. The 
later arrivals are associated with rays having larger launch 


and arrival angles and thus are subject to interaction with 


the sea surface. 
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Figure 26 Relative amplitude of all rays from source (A) 
arriving at a vertical surface 50 km away at location (B) 
for a sea state of 3. 


Е. RESOLVABILITY ANALYSIS 
1. Introduction 

Two methods for resolving ray arrivals using a 
vertical array of 16 hydrophones with a 10 m interelement 
spacing were considered. The first method used the 
hydrophones as independent omnidirectional receivers with 
resolvability in time. The second method used all the 
hydrophones together as a plane wave beamformer with 
resolvability not only in time but also in arrival angle. 
Thus, in order to estimate the total number of resolvable ray 


arrivals analyses of time and angle separations of the 
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eigenrays were required. For the time analysis the eigenrays 
for each individual phone were compared to each other for a 
minimum time separation requirement. For the angle analysis 
the plane wave beam pattern of the line array was used to 
construct the beamformed arrival structure. The beamformed 
arrivals within each beam were then compared to each other for 
the same minimum time separation requirement as used for 
individual phones. 

In а о to addressing resolvability it was also 
important to consider the accuracy of the travel time 
measurements. A resolved arrival may not be a useful datum 
for tomographic inversion if its travel time uncertainty o, is 
large. This uncertainty depends on the signal-to-noise ratio 
(SNR) associated with the individual arrival and can be 
estimated by considering the source level, ambient noise, and 
array gains from directivity and signal processing. Estimates 
of o, for both methods are included in this section. (Chiu, 
1992) 

2. Method I - Individual Hydrophones 

a. Time Analysis 
The investigation of resolvability in time was 
accomplished by comparing each ray arrival in time with its 
neighbor. If the separation in time was more than the width 
of the transmitted pulse then the ray was resolved from all 


other arrivals. The 224 Hz source signal has a bandwidth of 
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16 Hz. Therefore, the theoretical width of the received pulses 
mommy l6 HZ or 62.5 msec. 
b. σι Analysis 
Dnoxdcteminectheog "s of the resolvable arrivals an 
SNR estimate was needed for each arrival. Using an approach 
similar to Spindel (1979) and from Miller (1992) the SNR 


estimate at the receiver for a single pulse used was: 


SNR = RAMP - DNL + PCG + CAG (9) 


where RAMP is the ray amplitude (source level (SL) minus the 
total transmission loss (TL,,)), DNL is the detected ambient 
noise level, PCG is the pulse compression gain, and CAG is the 
coherent average gain. The SL for the 224 Hz source is 183 


dB. TL,, was computed for each ray as described in section E. 


t 
DNL was computed from the following equation (Kinsler et al., 


1982): 


DNL = NSL + 10log(BW) - DI (10) 


where NSL is the ambient S level, BW is the band width of 
the transmitted signal and DI the array directivity index. 
NSL has many contributions including shipping, agitation of 
Esa surface, bioacoustics, sea cm. and seismic sources. 
Ambient noise due to shipping and the sea surface was only 
considered here. For medium shipping density and a sea state 
of 3, a NSL of 67 dB was found in Kinsler et zu. (1982). The 


ambient noise information from Kinsler et al. (1982) is for. 
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deep water and was used here as an approximation. Using a 16 
Hz BW the DNL increased by 12 dB. For individual 
omnidirectional hydrophones DI = 0 so DNL = 67 + 12 = 79 dB. 

The PCG and CAG are signal processing gains. Using 
equations from Spindel (1979) and the planned signal 
parameters for the 224 Hz source, the following values were 
calculated. With a kernel sequence of 63 digits the PCG is 
10109(63) ог 18 dB. With 30 repetitions of the sequences, a 
CAG of 10l0g(30) or 15 dB was obtained. 

According to Spindel (1979) a reception with a SNR of 
about 20 dB is desired in order to clearly distinguish 
multipath arrivals from noise and to reduce the erro g 
travel time estimation to a figure adequate for tomographic 
inversion. Since the error is approximately given by 


(Spindel, 1979) 


ΤΙ. ---. (11) 


a 20 dB SNR and a pulse width of d = 62.5 ms yields a oot 
6.3 ms. Figure 27 provides a plot of SNR verses o. for Έσω 
values of d. 

Solving Equation (9) for RAMP with SNR 2 20 dB, DNL - 79 
dB, PCG = 18 dB and CAG = 15 dB yielded a RAMP » 66 dB. The 
implication was that resolvable arrivals with RAMP > 66 dB 


have travel time uncertainty < 6.3 ms. These strong arrivals 
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Figure 27 SNR versus o, for two pulse widths. The 224 Hz 
source has a d of 62.5 ms whereas the 400 Hz source has a d 
of 10 ms. 





would give high quality travel time data for tomographic 
inversion. 
c. Resolvable Rays 

To assist in identifying the resolvable гау 
arrivals a table summarizing all of the eigenrays for each of 
the sixteen hydrophones was compiled using a MATLAB (The 
Mathworks, Inc, 1989) computer program 'resanl.m' (Appendix 
A). The table is given in Appendix B and the list is in order 
of increasing travel time. For all the rays traced to a range 
of 50 km, the rays that arrived within 1/2 wavelength of each 


hydrophone were picked as eigenrays. 
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Each arrival was then compared to its adjacent arrival in 
time. Arrivals with greater than 62.5 ms separation to the 
arrivals before and after were considered resolvable in time 
and are annotated with a 'R' next to the arrival time in the 
table in Appendix B. 

The arrivals resolvable in time were then inspected for 
adequate g,. Arrivals with a o, < 6.3 ms corresponding to a 
RAMP 2 66 dB were considered as having adequate travel time 
accuracy (see Sec. F.2.b). Inspection of the table in 
Appendix B yielded 42 resolvable arrivals that met this 
criterion. Table III provides a summary of these 42 arrivals. 

Figure 28 shows the simulated arrival structure at the 
eighth hydrophone from the top of the vertical hydrophone 
array. This phone is located at a depth ot 220 m- The piot 
was created using the MATLAB program 'hyd8beam.m' (Appendix 
A). The figure shows a total of 27 ray arrivals at this 
hydrophone. The arrivals resolvable in time and have a o, < 
6.3 ms are annotated by R17, R18, R19, and R20. The numbers 


correspond to the numbering in Table III. 
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Table III INDIVIDUAL HYDROPHONE RESOLVABLE RAY ARRIVALS 


WITH o, s 6.3 ms 

Launch Arrival Arrival 

Angle Angle Hyd RAMP с, Time 

No.  (deq) (deq) No. (dB) (ms ) (s) 
1 11.49 -10.7803 1 82.0474 1.0 34.8858 
2 12.74 +10.6750 1 74.9298 2.2 34.9759 
3 14.83 +13.4689 1 72.2280 35] 34.3221 
4 11.51 -11.0427 2 80.6225 Lo 34.8846 
5 12765 +10.8853 2 82.0131 1 0 34.9881 
6 7.78 +8.0024 3 70.2431 3.9 34.7172 
7 1152 -11.1197 j 76.0964 2:50 34.8853 
8 14.74 +13.7346 3 76.0389 2.0 9009332 
9 8.62 + 8.2148 4 85.9763 0.6 34.7075 
10 11553 =11.2162 4 78.5257 5 34.8852 
11 14.60 +13.9993 5 77.0836 1.8 3503518 
T? 11.54 -11.3108 6 75.0689 2.2 34.8864 
13 14.53 +13.9967 6 76.7759 1.8 35. 3615 
14 11.55 -11.5583 7 76.8271 1.8 34.8868 
15 14.45 +14.0910 7 76.9178 1.8 35.3719 
16 13.10 -14.2447 7 66.9958 5.6 39.9378 
17 8.50 +8.4845 8 1925255 Пе 34.7140 
18 11.63 -11.3509 8 75.6184 221 5 06011 
19 14.38 +14.0919 8 77.3569 y 35.3814 
20 17.01 +16.5614 8 66.5519 5.9 35.8676 
21 7.68 -8.3897 9 ΤΘ. 1157 1.6 34.6960 
22 11.56 -11.3548 9 75.7405 20 34.8880 
23 11.61 -11.3504 9 73.5302 2196 35.0963 
24 14.29 +14.1118 9 78.1330 16 35.3934 
25 16.92 11ο 55955 9 67.5459 S2 35.8819 
26 14.18 S14 1081 10 912232 1.4 35.4087 
27 11:13 -14.2383 10 7372166 Zt 235.5376 
28 16.83 τις 5507 10 66. 6980 5.8 35.8963 
29 14.06 +14.0723 Ll 78.8534 1.4 35.4251 
30 13:15 -14.1958 1 2707723 2:0 35.5364 
31 ΠΟΘ πο ο η ΤΙ 66.8359 Տ» 35.9145 
32 Τα] -8.3325 12 0093256 1.4 34.6980 
ES 13.94 *14.0684 T2 78.8688 1.4 35.4415 
34 13:718 -14.1533 12 η 5155 29 35.5344 
BS 16.65 “57 12 66.1088 6.2 3579256 
36 172 -8.3059 13 83.3586 0.9 34.6982 
37 155-9 3 *14.0650 το 78.6391 υ 5 35.4564 
38 13.19 -14.1535 13 70.6539 3157 B5. 5337 
39 16.57 416.5326 13 66.5419 5.9 35.9380 
40 7.74 -8.2335 15 80.9442 Ἰώ] 34.6987 
41 12510 -11.2190 15 72.4460 3.0 35:0256 
42 8.45 -8.3734 16 74.1797 2.4 34.7021 
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Figure 28 Simulated arrival structure at hydrophone #8 
(depth 220 m) at a sea state of 3. 


3. Method II - Plane Wave Beamformer 

If the vertical array of 16 hydrophones is used as.a 
plane wave beamformer then angle, in addition to time, becomes 
another parameter available for resolving arrivals. To 
determine the resolution in angle, the plane wave beamforming 
method of Ziomek (1985) for a linear array of equally spaced 
point hydrophones was used. The far-field directivity 
function or beam pattern of the linear array can be computed 


uSing a Fast Fourier Transform (FFT) computer algorithm. 
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Using a MATLAB computer program ‘beampat.m' (Appendix A) 
incorporating such FFT algorithms, the theoretical beam 
pattern of the hydrophone array was computed. 

Figure 29 shows the beam pattern of the vertical receiver 
array with the beam steered to 20° off broadside. The 
corresponding beam width defined by the 3 dB down point 
(0.7071 normalized directivity) is about 3.0°. The beam 
pattern steered to broadside was found to be about 2.5° wide. 
The grating lobe for this array appears 48° away from the 
steered beam so should not cause a problem for arrivals within 


20° of broadside. 
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Figure 29 Beam pattern for modeled linear array with beam 
L eered to 20 degrees off broadside. 


To assist in identifying the resolvable beam formed 


arrivals, a table was compiled using 'resanl.m'. This table is 
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presented in Appendix C. The method used here was similar 
Method I for individual hydrophones with the following 
exceptions. Following extraction of all eigenrays arriving 
within 1/2 wavelength of each hydrophone, the eigenrays were 
separated into fifteen groups associated with their arrival 
angles. The 15 groups are associated with 15 nonoverlapping 
angular sectors. Each angular sector is three degrees wide. 
The chree degree width was based on the beam pattern analysis 
discussed above. The limiting arrival angles defining these 
sectors are given at the end of the table in Appendix C. 

Similarly, the beamformed ray arrivals were compared in 
time for a 62.5 ms separation. The rays meeting this time 
separation criterion within each angular sector are annotated 
with an 'R' next to the arrival time in the table in Appendix 
er 

For the c, analysis, the same parameter values foP PCOR 
CAG and BW as before are used here. The only difference is 
that a directivity gain was added. The directivity index 
(DI), according to Ziomek (1985), may be approximated . by 
lOLog(number of hydrophones). For 16 phones it is 12 dB. 
Thus from Equation (10), DNL is now 67 + 12 - 12 = 67 GB. 
Using Equation (9), this results in a minimum RAMP requirement 
of 54 dB for resolvable arrivals 1 ατα. πα]ΠΕν» 5Η 
6.3 т. 

By inspection of the table in Appendix C, 11 resolvable 


beamformed arrivals that also have an adequate c, of 6.3 ms or 
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less (i.e., with RAMP 2 54 dB) were obtained. Table IV 
provides a summary of these 11 arrivals. Seven of these 


arrivals are unique to the 42 arrivals determined using Method 





T. 
Table IV PLANE WAVE BEAMFORMED RESOLVABLE ARRIVALS WITH σι 
< 6.3 ms 
Launch Arrival Arrival 
Angle Angle RAMP σι Тїте 
Мо. (дед) (deq) Sector (dB) (msec) (sec) 
1 15.59 -16.5311 2 62.4927 2.4 36.0685 
2 15.41 -14.2058 3 55.6390 5, 2 35.2809 
3 13.10 -14.2447 3 64.9958 Ig 3525378 
4 11.56 -11.3548 4 74.7405 0.6 34.8880 
5 13.02 -12.2672 4 733231 0.7 35.0269 
6 ПБЯ -11.3504 4 73.5302 0.7 35.0963 
7 1171 +11.1330 12 67.3867 1.3 35.1003 
8 14.83 +13.4689 12 72.2280 0.8 359592299 
9 13.55 +14.0571 13 77.6399 0.4 35.4936 
10 29 +16.3884 13 64.2685 1-9 358375 
11 16.35 +16%”%5172 14 65.0820 178 35.9726 
i 
Figure 30 shows a total of three arrivals in angular 
sector 13. Arrivals resolvable in time within this sector 


which have a DENS 6.3 ms are annotated with 'R9' and 'R10'. 


The numbers correspond to the numbers in Table IV. 
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Amplitude (Pa) 


35.7 35.8 
Arrival Time (sec) 


Amplitude (dB) 


ա», JEUX 
Arrival Time (sec) 





Figure 30 Simulated beamformed arrival structure in sector 


13 (+13.5° grazing angle). R9 and R10 are resolvable and 
have o, < 6.3 ms. 
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VI. CONCLUSIONS 


The objective of this thesis was to examine the arrival 
structure of the planned tomographic transmissions in the 
upcoming BSPFEX using computer simulation. This study 
examined the resolvability of those rays that intersected a 
vertical hydrophone array 50 km from a near bottom source. 
Conclusions are: 

1. The duration of the multipath arrival structure is about 
4 s. 
2. The arrival structure can be grouped into the followinq 


two distinctive regimes: 
a. Early Arrivals (0-600 ms) 

(i) These arrivals have angles less than 10?. They 
exhibit the least transmission loss due to minimum 
interactions with the sea surface. 

(11) They exhibit mode-like properties. Large 
‘number of simultaneous ray arrivals аге found. The 
investigation of the resolvability of modes will require a 
simulation using full-wave acoustic models based on parabolic 
equation or coupled mode methods. 

(iii) Fewer ray arrivals are resolved in this 


regime. 


Ээ 


b. Late Arrivals (600-4000 ms) 
(i) These arrivals have angles greater than 10°. 
They exhibit higher transmission loss due to greater 
interactions with the surface. For sea states 25 (i.e., mean 
wave heights > 4 m) little signal amplitude is expected of 
rays launched at > 13°. 


(11) Individual ray arrivals are separated more in 


time. 
(111) More ray arrivals are resolved in this 
regime. 
3. One to four accurate ray arrivals are resolvable per 
Single hydrophone. Accuracy is measured by o, and a 


requirement of o, s 6.3 ms is used for extraction of the 
accurate arrivals. 

4. A total of 42 ray arrivals with the required amplitude 
are resolvable from the sixteen individual hydrophones. All 
of these arrivals come within the first 1500 ms in the 
multipath arrival Structure. 

5. Seven additional, non-redundant, ray arrivals with the 
required amplitude are resolvable through plane wave 
beamforming using the 16 hydrophones. All of these arrive 


within two 6 degree bands of +11 to +17° and -11 to -17?. 
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APPENDIX A: LIST OF EXTERNAL PROGRAMS 


The following programs were developed or extensively 
modified during the course of this research. They include: 
ւ. cordat.f - Path and time correction for a vertical 


receiver and TL calculation. 


2. tgridin.f - HARPO topography input generator. 

9. gridin.f - HARPO sound speed profile input generator. 

4. beampat.m - vertical hydrophone array beam pattern 
generator. 

5. resanl.m - Ray resolvability analysis table generator. 

6. hyd8beam.m - Arrival structure plotting program 


including pulse shape generation. 


55 


Jun 8 08:37 1992 cordat.f Page 1 


С cordat £ 
С 


σςσσσσσσσσσσσσσσσσσσσσσσσσοσσσσσσοσοσσσσσσσσσσσσσσἒσσσσἒσσἒσσςσἒσσσσσςσσςἒς:,,. ΜΕ 

С 

с Written by: J. Mark Elliott, LCDR/USN, NPS, Dec 1991 

C 

ссссссссссссссссссссссссссссссссссссссссоссссссссессссссссссссессессс ыыы ше ыыы 
Modified 23 Jan 92 by J.M. Elliott to correct for total ray travel time 

to the vertical receiver surface. 


Modified 27 Jan 92 by J.M. Elliott to calculate Transmission Loss due to 
spreading. 


Modified 21 Feb 92 by J.M. Elliott account for trans loss due to bot and 
surface reflections. 


Modified 20 Apr 92 by J.M. Elliott to count the no of surface, bottom, 
and turning point interactions. 


ՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇԵՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇՇԵՆՆԵ ՆԵՆ ՆԵՆԵՆ 


ռո գ ооо о ооо 


Purpose: 

This program extacts from RAYSET (HARPO output file) the last line of 
data for a particular launch ray corresponding to the first event (ray point) 
after passing max range (W28 of DINP). Then this is backed up the ray path 
geometrically to compute the height above msl and arrival time for a ray 
crossing a vertical recvr surface at max r. 

This program also computes the fol: Transmission Loss due to spreading, 
surface interactions, bottom interactions; and counts the number of surface 
interactions, bottom interactions, and turning points. 

This program will skip rays that are trapped above the surface. This 
problem was later found to be mitigatied but not eliminated by reducing W44 
(initial integration step size) from .lkm to .05km or less. 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


Define the 10 Columns of RAYSET data output, max range specified in 
DINP W28 (rmax), and initial launch elevation angle specified in DINP 
W15 (ths), and the step in launch ele angle Wl1?7 as dths. 

Also six indexes are used i, tlbot, sloss, botcnt, srfcnt, and tpont. 
The index i is to be used in the spreading trans loss calc. The indexs 
tlbot and sloss are used to measure the TL due to bottom and surface 
reflections respectively that the ray makes along its total path. 

The number of bottom, surface and turning points are'counted by 'botcnt', 
'srfcnt',and 'tpcont' respectively. 


O000000000000000000000000 


real*4 £f(lÜü),ths,thr,r,zp,duml,t,tlbot,sloss ЕТЕ 2 rp ds ο ο 


a cavel,cave2,z2,rmax,rmax2,thr2,h,all,a1l2,a21,a22,b1,b2, 
a w, zl,rl,tlspn,tlspo,theta,sigma,k,botcnt,srfcnt,tpcnt 
rmax = 50. 
ths = 0.00 
dths = 0.01 
i = 1 
tlbot = 0 


sloss = 0 
botcnt = 0 
sri cnt 0 
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tpcnt = 0 
Ë Define output files: 
[с raysetX Y - rayset data from harpo output, launch angles from X to Y 
C Ε ταας - launch angle vs arrival depth output 
ο f3time.dat= arrival time vs arrival depth output 
С f3tl.dat = launch angle vs transmission loss due to spreading 
Ic f3ttl.dat = arrival time, total trans loss (dB) 
€ mres dat = ths, thr, z, tpent, srfcnt, botcnt, ramp, t 


Open (unit=10, file=’ rayset0_10’ ,status="°old’, form 


A =' formatted’ ) 
Е open(unit-20,file-'f3.dat',status-'new',form-'formatted') 
C open (unit=30, file=’ f3time.dat’,status=’new’, form=’ formatted’ ) 
С open(unit-40,file-'f3tl.dat',status-'new',form-'formatted') 


open (unit=50,file=’f3ttl.dat’, status=’ new’ , form=’ formatted’ ) 
open (unit=60, file=’ f3res.dat’,status=’ new’ , form=’ formatted’ ) 


C skip the first 7 lines of RAYSET 
PO I5 i - 1,7 


read(10,9*) 
15 continue | 


Read in the 10 columns of RAYSET data and call the data n for new 
Define RAYSET columns: 
#(1) = earth radius to mean sea level (msl) in km 
£(2) = earth radius to ocean floor in km 
#(3) = horiz range to ray point along msl in km 
Ε(4) = earth radius to ray point in km 
£(5) = ray phase angle 2 
f(6) = ? 
f(7) = phase arrival time in sec 
f (8) = local sound speed in m/s 
£(9) = 3 
£(10)= local ray elevation angle in deg 
1 ΙΙ. ЕСЕ Է(4),6 (5) է (6), 207). 208), (9) , Е (10) 
Еп = ք(1) 
f2n = £(2) 
fom = £ (3) 
ք4ո = £(4) 
Hn = £ (5) 
աջը = £ (6) 
En = Е(7) 
fen = £ (8) 
f9n = £(9) 
f£lOn = #(10) 


Ω 


C 
С Look for the last line -999 flag of each launch ray data set, otherwise 
C (else) call the new data line the old (fxo) data, compute surface 
c and bottom trans loss, and then look at the next data line in RAYSET: 
ο πι τε 900. )then 


> 
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go tous 
end if 


Compute bottom/surface interaction trans loss. Let each bottom ref 
equal 0.5 dB of loss. For surface loss compute the trans loss 
(sloss) using gaussian scattering theory (Clay & Medwin, 1985): 
Define: 

Zp = arrival ht above msl in km 

duml = arrival ht above bottom in km 

sloss = ~20*1l0g10 (abs (exp (-2* (k**2) * (sigma**2) * (sin(£(10) ) ИЕ 

k = wavenumber at surface using avg surf sound speed of c(m/s) and 

frequency f (Hz) 
- 2*pi*f/c -2*pi*2247146S RON SD 

sigma = rms wave ht in meters( 0.7071 used for sea state 3) 
0.0707 for sea state 1 
0.7071 for sea state 3 
2.8284 for sea state 5 
graze angle assumed to be appx = local ele angle f10 
botcnt <= no of bottom interactions 
srfcnt = no of surface interactions 
tpcnt “по of turnimg Ponts 


l Ww H l 


G) GO G G о Ооо осо ооо ооо ο ο GO оу оо 


First test to see if the RAYSET data line was a repeat by checking range £3, 
and if so skip the computation and repeat (cause of this problem is unk): 
1f (f3n.eq.f3o)then 
ао Ες 1 
end if 
zD m= гап = CIN 
duml - f4n - £2n 
if(abs(dum1).1e.0.00009) then 
Έ]δος -]σσος ο 
botcnt ου πο. 
end if j 
С k = 0,960; 
С Sigma = 2.8284 
if(abs(zp).le.0.00009) then 
theta = sin(fl0n*3. 142597 130) 4-2 
С sloss = -20*1oglÜ0(abs(exp(-2*(k**2)*(sigma**2)*Eheta) DONNE NC 
c Can't get program to compute sloss in this form, cuse unk, must осеке 
Sloss = -20*log10(abs(exp(-0.9229*theta))) * sloss 
srtent - srtentes e 
end if | 
Now count the no of turning points. Since the method I use just counts the | 
no of times grazing angle changes sign I subtract out the sign changes due to. 
surface and bot bounces at the end of this program just before the write 
statement: 
if (abs (£100-£10n) .ge.abs (£100) ) Ehen 
ἵροου = Cpcn tat. 


Q © 


о Оооо 


end if 
С 
Փիր 
f20 = Էշը 
f3o = f3n 
f4o = f4n 
ԷԹ» fon 
f60 = fon 
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ETO ETN 
£80 f8n 
ք9օ f9n 
ք109 » fl0n 


Soto 1 


Look when f(3)(range) » max range specified(W28) otherwise (else) the ray 
is trapped above the surface and must be skipped 


o o o ԹԹ 


2 if(f3o.GT.rmax)then 
ШО to 3 
else 
Шис = ths + dths 
go to 1 
end if 


Correct the ray arrival depth back up the raypath from the last event 
to where it crossed an imaginary vertical surface at r max (W28 of DINP). 
զ -Ջ ուս, ԱՔՏ, 1991). 

Define 
arrival ht above msl (km) 
arrival ht above bot (km) 
local arrival angle(deg) 
msl dist from source(km) 
msl dist to recvr (km) 
ray arrival hieght above msl accross a vert sfc at 
the recvr (km) 

ths = launch angle last + step in elevation angle(W17 DINP) 

3 Zp = f4o - flo 


Zp 
duml 
Enr 
rp 

թ 

z 


QOO O ОООО О ՈՌ СОГ ООО 


duml = f4o - f2o 


pur - flOo 
το - Εσο 
r = rmax 


c If the ray has just reflected from the bottom or surface then correct the 
Beomagte to incident vice reflected by changing the sign: 
if (abs(zp).le.0.00009.0r.abs(dum1).1e.0.00009) then 
enr -thr 
end if a 
HL: r = rp) * tand(thr) 


e արատ 20.100) Էհտ, 2 

© ШЕ ο πας (2Ε9. 4) 

С ο τν (for ray travel time: 

e Define: 

C ds = path length past vertical receiver surface 

Ç Ер = total phase travel time to horz recvr surface 

e Git — travel time for path beyond vert recvr surfc 

թ fe = corrected total ray travel time to vert recvr surf 
С cavel= ave sound speed for 0-150m array depth = 1.460km/s 
С cave2= ave sound speed for 150-320m array depth=1.446km/s 
С 


cavel = 1.460 
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cave2 = 1.446 


tp = f7o 
ds * sqrt((z = zp) 72 Fr =- τοι o) 
dt, = ds/cave2 


t - tp - dt 
wrritet30 200072 πε 
200 format (2£9.4) 





Calculate TL due to spreading (Ching-Sang Chiu, NPS, 1991): 
Define: 
ths launch angle 
ths2 previous ray launch angle | 
E dths = difference btn one ray launch angle and the previous one 
2 arrival depth 
22 previous ray arrival depth 
rmax vertical reciever range 
rmax2= rmax (previous vert recvr range) 
ЕВЕ arrival angle 
ЕВЕ 2 previous ray arrival angle 
h xsec distance btn launch ray and previous ray, perp to 
launch ray at rmax 
tlsp = trans loss due to spreading (dB) 


QQ Oo Q QQ O O O O 0 0 00 090 0 0 0. 0. 0-0 


it uedg σου ας 5 


Ω 


if (abs(thr).1e.0.00009)then 
h = abs (z-z2) 

go S 

endif 


if (thr.gt.0.00009) then 
phi = -( 90. - &һг) 
endif 


if(thr.1t.-0.00009)then 
phe + Icha 
endif 


all = 1. 

а12 -tand(thr2) 

а21 РАИ 

а22 -tand(phi) 

z2 - rmax*tand(thr2) 

Ե2 z - rmax*tand(phi) 

ч = all*a22 - a21 * a12 

zi  (bl*a22 - b2 * a12)/w 

τι “πιο: 

h = sqrt((z =) Zh) **~2 2) ү 2592 1-7 


| կ 


у 
к— 
կ | 


3 


c Check for Caustic by looking for h near zero. TE Տօ ՅՏՏւզոը the rr Шш 
c loss cale from the last ray (tispo) еда ԵԱ Երի came Έντα. 
if (abs (h) .le.0.00009) then 
tlspn = tlspo 
GQ του” 
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end if 
6 ա բո = 10*10g10( (rmax*h*1000**2) / (cosd(ths) *dths*3.14159/180.)) 
5 22 - 2 

mice = thr 

ths2 = ths 

rmax2 = rmax 


write (40,300) ths, tlspn 
300 format (2f£9.4) 


Compute total trans loss (tlt) for each ray as sum of bottom loss 
(tlbot), surface loss (Sloss), and spreading loss (tlsp). Trans loss 
due to absorption will be neglected since its on order 0.3dB for 224 Hz 
over 50 km distance: 

ШЕ = tibot + sioss + tlspn 

Compute relative amplitude (ramp) by subtacting the source level (183 dB) 

from the total transmission loss (tlt): 

ramp = 183 - tlt 


Ее (50.400) ths, thr, z, t, tit, tlbot, sloss, tlspn 
400 format (8f£9.4) 


ա ԱՆ = tpcnt - (srfcnt + botcnt + 1) 
ff (tpent .1t.0.00009) then 
tpent = 0 
end if | _ 
ԱԱ 50500) ths, thr, z, tpent, srfcnt, botcnt, ramp, t 
ա աթ Հ Էոմե (889.4) 


i = i + 1 

ths = ths + dths 
tlbot - 0 

Sloss = 0 

tlspo = tlispn 
tpcnt = 0 

srfcnt = 0 
botcnt = 0 


Repeat the process for the next launch ray ‘data set 


ао το 1 
close (10 
close (20 
close (30 
close (40 
close (50) 
close (60) 
end 


) 
) 
) 
) 
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Written By: John M. Emblidge Lt/USN, NPS, May 1991 
Purpose To write the manually gridded topography of the Barents 
Sea Acoustic Tomography Transmission Test study area to 
to an output file (GRIDBATHY.DAT) for use by tgridder.f 


bathymetry area closer to Bear Is in support of the 
Barents Sea Polar Front experiment to be conducted Aug 92 


J.M. Elliott, LCDR/USN, 16Jan92, 


grid to a true N-S coord sys. 
k * k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k hehe oe ode oe oe oe oe oe oe ode oe oe oe oe oe oe oe oe oe ooo A £ k k k <x x 


* 
k 
* 
* 
* 
* 
Modified by: John M. Elliott, LCDR/USN, 28DEC91, to change location ΘΕ. 
* 
* 
* 
to correct the Dachy k 
* 
* 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 


INTEGER I,J 
LOGICAL EX 


PARAMETER (MX=13, MY=13) 

REAL*8 BLAT, TLAT, DELLAT, DELLON, LLON, RLON 
DIMENSION BOTGRID (MX, MY) 

DIMENSION XCOORD (MX) , YCOORD (MY) 


DIMENSION 111 (47), 182 (MY), IN3 (MY) , IN4 (MY) , INS (MY) , IN6 (MY), 
! IN/ (MY) , INS (MY) , INS (MY) , INLO (MY) INT I (MY P INI ee 


IN13 (MY) 
e ЇМ14(МҮ),ТМ15(МҮ),ТМ1б6(МҮ),ТМ17(МҮ),ТМ18(МҮ), 


© ЇМ19(МҮ),ТМ20(МҮ),ТМ21(МҮ) 


k k k kk k k k k k k k k k k k k k k k k k k k k k k k k k eode obe ode che oh oe ode oe eoe oe ode che oe oe ode oho ode oe oe oe oe oe oe ode oe oe eoe ode ko e ook k wk * 


C DEFINE BATHYMETRY DAT. TOP ROW IS EASTERN BOUNDARY, 
C VERTICAL ROW IS NORTHERN BOUNDARY: 


RIGHT HAND 


DATA IN1 /219,220,175,137,119, 116113 το Е 
DATA IN2 /310,245,201, 159,130, 118111, 121 12121277 նն Ա 
DATA IN3 /340,210,226,190,156,139,130,1352 12212211 ՄԱ ն 
DATA IN4 /359,295,250,222,192, 169, 152,145) e776 --. 
DATA INS /369,316,275,249,220,194,173,160,150, 111. 115 ΤΙ... 
DATA IN6 /386, 340,302,272,241,210,197, 184, 166,156,153) KDP 
DATA IN? 4/406, 360,323,296, 260,224, 210,203, 182, 167, 1567 ioe 
DATA INS8 /420,382,345,312,219,243,220201 19218151 ն 
DATA IN9 /428,401, 366,330, 298,268, 240, 217, 2097 ο 16 P l s DD 
DATA IN10 7/438,411,375,338,312,2907265/23987 220 478 0 NT) n O 1 
DATA IN11 /442,419,384,342,322, 306,288,266, 2407227722372 a ee 
DATA ԼԱ112 /440,422,393,3583312316 03222222 22. Ն` 0 
DATA IN13 /433,421,400310,340,325315 50242012: 9 Ե 


X Ck ce ck oe eee eoe oe ce eee e ce eoe oe eoe eoe oe eoe e eoe oe e oe eoe oe oe oe eoe oe oe eoe e ce oe eoe ce oe oe eoe oe oe oe oe eoe dee ove x x kx kx x X 


C DEFINE THE GEOGRAPHICAL LIMITS OF BATHYMETRY BOX: 
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Pee LAT, TLAT, LLON,RLON / 74.22d0,74.97d0, 22.72d0,25.52d0 / 


` k k k k k k k k k k k k k k £ k k k k k k k k k k k k k k k kk kk kk kk xx xk xk k kk kk kk kk kk kk kk k'ik kk kk kk k xk kk kk kk k 






2 Շ TO SEE IF THE OUTPUT FILE ALREADY EXISTS, AND IF SO DELETING 
> IT AND CREATING A NEW ONE TO WRITE THIS RUNS DATA TO. 


INQUIRE (FILE-'GRIDBATHY.DAT',EXIST-EX) 

IF (EX) THEN 
OPEN(35,FILE-'GRIDBATHY.DAT',STATUS-'OLD' , FORM=’ UNFORMATTED’ ) 
OPEN (35,F ILE=’ GRIDBATHY .DAT’ , STATUS=’ DELETE’ , FORM=’ UNFORMATTED’ ) 

ENDIF 

OPEN (35, FILE=’ GRIDBATHY .DAT’ , STATUS=’ NEW’ , FORM=’ UNFORMATTED’ ) 

REWIND 35 


k k k te He He He KK HK HK KKK KK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KK 


© CALCULATING THE DEGEE SPACINGS BETWEEN LAT AND LONG GRID LINES. 


He He HK KK KKH KH KH KH KK KKK KKK HK KK KK KKK KK HK HK KKK KKK KKH KKK HK KK KEKE KKH KEKK KKK ko 


DELLAT 
DELLON 


(ТАТ - BLAT) / (MY-1) 
(RLON - LLON)/ (MX-1) 


ELFILLING THE X & Y COORDINATE ARRAYS 


XCOORD (I+1)=LLON -* DELLON*dfloat(I) 
ENDDO 


DO J=0, (MY-1) 


YCOORD (J+1) = BLAT + DELLAT*dfloat (J) 
ENDDO 


k k k k k kk k k kk Kk K Kx kx yk ki k kk kk xk xk kk xk eoe ok ok eoe oe oe dece oe oe ok ok kk k kk kk kk kk kk kk kk kk kk kk kk 


PUPFORMING THE BOTGRID MATRIX TO BE WRITTEN TO THE OUTPUT FILE. 


DO I21,MX 
DO J-1,MY 
Ir I EO. 1) BOTGRID(I,J) = -ԼԱԼ(Մ) 
IF(I .ЕО. 2) —BOTGRID(I,J) *  -IN2(J) 
աթ Ե ο] ΕτοΕτω == —IN3(J) 
πο ο ο 4) BOTORID(I,J) = -IN4 (J) 
ει ο ο) BOTGRID(I,J) - =INS (J) 
προ ο ο) SOTGRID(I,J) = —-IN6(J) 
BREMEN E DOTORID(I,J) =  -INT7(J) 
ποτ εἰ Ὁ ΡΟΤΟΕΙ ΡΕ ΤΠ) = -IN8 (J) 
ПЕ ᾿ν, ΕΟΟ J) = -—IN9(J) 
το ο ο. 1ο) BOTGRID (IMT) = —INIÓO(J) 
EXON ОЗЕ ВОТСЕТО (217) = -IN11 (J) 
πι ο. աշ) BORTCRID(I, J) =  -IN12(J) 
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IF(I .EQ. 13) BOTORID( (T Թ-ի 
ENDDO 
ENDDO 


k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 


C WRITTING XCOORD, YCOORD, & BOTGRID TO THE OUTPUT FILE GRIDBATHY.DAT 


WRITE(35) XCOORD 
WRITE(35) YCOORD 
Wrlte(6,*) xcoord 
write(6,*) ycoord | 
DO I = 1,MX 

DO J = 1,MY 

WRITE (35) BOTGRID (I,J) 
С write(6,*) BOTGRID(I,J) 

ENDDO 

ENDDO 





k k k k k k k kk k k k k k k k %@ k k k k k k k k f k k k k k k k k k k k k k k k k k £ k k k k k k k k k k k k k k k k k k k k k k k k k k k k k + 


CLOSE (35) 
STOE 
END 
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C gridin:f 


QU ec k k k k k k k k k k k k k k k k k k kk k k kk k k k k kk k e he e kk k e e e e e e ke e e e e e e e e e e v v ke ke ke ke ke e e ek 


Written by: John M. Emblidge, LT/USN, NPS, May 91 
FOr : Barents Sea Acoustic Tomography Experiment. 


Purpose : To create the input file (GRIDSVP.DAT) for the HARPO 
subroutine gridder.f. This file must include the 
spacings in the 3 othonormal planes, and the SSP at 
each point. 

The file is written so that gridder.f reads the data as 
follows: 

DINE Ts XCOORD (MX) 

LINE2: YCOORD(MY) 

LINE3: ZG(MZ) 

LINE4: SVP (MX,MY,M2) 


Structure : The most important thing to remember about using this 
routine is that the values of MX,MY,MZ must be the same 
as they are in GRIDDER, otherwise the data will not be 
read by GRIDDER as you intended it to be. 


Modified by: J. Mark ELLIOTT, LCDR/USN, NPS 
16 JAN 92: TO CORRECT FOR NEW BATHY in support 
of the Barents Sea Polar Front exp - 50km box 
centered at 74.6N 24.12E. 
20 Jan 92: to correct for new SSP field 
Orientation to along diagonals vice horizontals 


KKK KKK KKK KKK KK KK KKK KKK KKK KKK KKK KKK KK KKK KK KKK KKK KKK KKK KK KK KK KKK KK 


+ Ի + + + + HF F «2 Հ Հ Հ «ծ «էէ է Հ F OF F OF OF OF OF X 


DATA DICTIONARY : Definitions of each variable used 
CHARACTERS: 
INTEGERS 
REAL : BLAT = LATITUDE OF BOTTOM. OF BOX (SOUTH EDGE) 5 


DELLAT - NUMBER OF DEGREES BETWEEN GRID SPACINGS* 
DELLON = NUMBER OF DEGREES BETWEEN GRID SPACINGS* 
LLON LONGITUDE TOF LEFT SIDE OFTBOX ы 
RLON DONGLIUDE OF RIGHT SIDE OF BOX 

TLAT = LATITUDE OF TOP OF BOX (NORTHERN EDGE) 


* 

* 

ARRAYS : C1 = SVP IN FIRST REGION OF BOX * 
C2 = SVP IN SECOND REGION OF BOX * 

СЗ = SVP IN THIRD REGION OF BOX * 

CSPD = OUTPUT FOR GRIDDER * 

* 

PARAMETERS: MX = + 
Ж 

* 

* 


ՍԱ ԱՆԻ թո ΙΙ ΠΕ EOQUATER 
| 


OF GRID SPACINGS TP РЕЛЕ MERIDEN 


| 
| 
OF VERTICAL POINTS 


= 
ՒՀ 
li 

ՀԷ ՀԻ ՀԻ 


ΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΟΩ 
iS 
δ3 
| 


kk k k kk kk k k k k k k Kk k k k k k k k k k k k k k k k k K k k k k k k k k k k k k k k k K k k k k k k k k k k k k k k k k k k k k k 


IME CICITIDOUBLE PRECISION (A-H,O-2) 
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INTEGER L, J; K IX TI nT 
LOGICAL EX 


PARAMETER (MX = 13, MY = 13, MZ = 151) 
REAL*8 BLAT, DELLAT, DELLON, LLON, RLON, TLAT 


DIMENSION C1 (MZ) ,C2 (MZ) , C3 (M2) 

DIMENSION AVE1 (MZ) , AVE2 (MZ) , AVE3 (MZ) , AVE4 (MZ) , AVE5 (MZ) 

DIMENSION AVE6 (MZ) , AVE7 (MZ) , AVE8 (MZ) , AVES (MZ) , AVE10 (MZ) , AVE11 (MZ) 
DIMENSION XCOORD (MX) , YCOORD (MY) , 2G (MZ) 

DIMENSION CSPD(MX,MY,M2) 


DATA LLON, RLON, TLAT, BLAT / 22.72d40,25.52d40,74.97d40,74.22d0 / 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkŘkkŘkkkkŘkkkkkkkkkkkkkkkkkkěkkkkkŘkkkkkkkkkkkkkkkkkěkkěkkkkxk 


INQUIRE (FILE-z'GRIDSVP.DAT',EXISTZEX) 


IF (EX) THEN 
OPEN (74,F ILE=’ GRIDSVP .DAT’ , STATUS=’ OLD’ , FORM=’ UNFORMATTED’ ) 
OPEN (74,FILE=’ GRIDSVP .DAT’ , STATUS=’ DELETE’ , FORM=’ UNFORMATTED’ ) 
ENDIF 


OPEN (74, FILE=’ GRIDSVP .DAT’ , STATUS=’ NEW’ , FORM-' UNFORMATTED') 


REWIND 74 


հա k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 


C OPENING THE FILES THAT HOLD THE SSP'S CREATED IN MATLAB 


Kama KKK KKK Kah KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KK KK KK KKK KKK aK 


OPEN (41, FILE=’ ssplf.dat’, STATUS=’ OLD’ , FORM=’ FORMATTED’ ) 
OPEN (42, FILE=’ ssp2f£.dat’, STATUS=’ OLD’ , FORM=’ FORMATTED’ ) 
OPEN (43, FILE=’ ssp3f.dat’ , STATUS=’ OLD’ , FORM=’ FORMATTED’ ) 
OPEN (44, FILE=’zg.dat’, STATUS=’ OLD’ , FORM=’ FORMATTED’ ) . 
REWIND 41 

REWIND 42 

REWIND 43 

REWIND 44 


* 
www kw Km KK KK KK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KK KKK 


C READING IN THE THREE BASIC SOUND SPEED PROFILES, AND DEPTH ARRAY 


f k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 


ВЕАЮ (4105) (СПИ τι ο) 
READ (42, *) ;(G2 (1) yp I 1 MZ) 
READ (43,7) “(C3 (tye i= no) 
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READ(44,*) (ZG(I), I-1,M2) 


| 


K k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k £ k k k k 


C CALCULATING THE DEGEE SPACINGS BETWEEN LAT AND LONG GRID LINES. 


k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k +* 


DELLAT 
DELLON 


(TLAT - BLAT)/ (MY-1) 
(RLON - LLON) / (MX-1) 


| de ck ck e ke ee ok FFE Fe II TF FEI KIKI FO II FI I TO IK k k k k k k k k K k K k k k k K K K k k k k k k k k k k k k k k kk k k k k k k k K K 


С CALCULATING THE AVERAGED PROFILES TO FILE SMOOTH THE OVERALL C FIELD 
С WITH THE FRONT 


k k K k k k k k k k k k k k k k k k k k k k k k k < k K K k k k Kk k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k + 


DO I = 1,MZ 





աի 11)  (Cl(I) + с2(ї))/2 

ОЕ 1) թ» (Ը1(1) + AVE2(I))/2 

AVE4(I) - (C2(I) * C3(1)0)/2 

AVE3(I) = (C2(I) Հ AVE4(I))/2 

AVES(I) = (C3(I) + AVE4(I))/2 
ENDDO 


I ck ke oe kc oe oe ke ke e he eoe e eee e ke e ke eoe oe e ke ke ee ee ee ՆԱՆԱՆ ՅԱՆԻՆ ԻՆ ՅՆԱ ՆԵՆ Նանա KK 
| 


;C CREATING A NEW SET OF AVERAGES FOR THE FIELD WITHOUT THE FRONT 


KKK KKK KKK KKK KK KK KKK KKK KK KKK KKK KK KKK KKK KKK KKK IKK KK KKK KKK KKK KKK KKK KKK KKK 


DO I = 1,M2 


SED (JT) = (οἱ (τὸ * C3(I))/2 
| AVE7(I) = (C1(I) + AVE6(I))/2 
AVE8(I) = (C1(I) + AVE7(I))/2 
AVE9 (I) = (AVE6(I) + C3(I))/2 
AVE1O(I) = (AVE9(I) + C3(I))/2 
AVE11 (I) = (AVE1LO(I) * C3(I))/2 
| ENDDO 
ծտ"տ"տ"' EL E x X X kc kk OC ko KK OK e Kok ek e eK AC XO Je CK e KO e Յո: 
| 


[C FILLING THE CSPD ARRAY WITH VALUES 
'C ALSO FILLING THE XCOORD, AND YCOORD ARRAYS 


pk ke kk e koe ke ke e eee e koe oe eoe e eode e de oe oe ode eoe ode eoe de oe ode oe oe oe e eode eoe oe de ke ode oe oe de koe ke e ck koe e e kk ck eee ek oko kk oko 


To DO IX = 1,MX 

[ο ΄ DO IY = 1,MY 

Ё DO IZ = 1,MZ 

‘a ար (նր ՆՈՂ. 4) ορ Արո ) = C] (12) 
Te ШОО τν ООЛО ЕХ. LT. 6) CSPD(IX,IY,IZ) = AVE8(1Z) 
IC րար որ ԱՆ TIX. LT. 7) CSPD(IX,IY,IZ) - AVE?7(IZ) 
С րար ա ՆԵ IX. LT. 8) CSPD(IX,IY,IZ) = AVE6 (IZ) 
IC πο "πο И т. 9) Շտթը (14, 17, 12) » "ՄԲ9(12) 
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© IF(IX. GT.8 .AND. IX. LT. 11) CSPD (tu IY, IZ) ο ΡΕ... 
e IF(IX. GT.10 .AND. IX. LT. 12) CSPDUIX,1Y,12Z) = Avalos 
e IF(IX .СТ. 11) ՇՏՔը (12, 171127 1Թ-:155 012) 
G write(6,*) CSPD(IX ТУХ) 
C ENDDO 
Е ENDDO 
C ENDDO 
DO I=0, (MX-1) 
XCOORD (I+1)=LLON + DELLON*d£f£loat (I) 
C write(6,*) XCOORD (I) 
ENDDO 
DO J = 0, (MY-1) 
YCOORD (J+1) = BLAT + DELLAT*dfloat (J) 
С write(6,*) YCOORD(J) 


ENDDO 


*k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k Kk 


C ASSIGN SSP'S TO DIAGONALS OF BATHYMETRY COORDINATES - ORIENTATION IS 
C ALONG BOTTOM LEFT TO TOP RIGHT DIAGONALS: 
k ck de e kc e kc ke kc e e de e e ke e e e ke ke oe e ce e e e ce ce ke ee e e ek e e e kc eoe e ke e ke koe e e kc ke ke e e oe ee e kc ke koe ee e e ke ke ek 
COUNTERY 5 2 
DO IX = 1, MX-2 
DO IY * 1, MY-COUNTERY 
DO IZ = 1, MZ 
CSBDIGEX, IY, I2)" :9630052) 
ENDDO 
ENDDO 
COUNTERY =“ COUNTER ta 
ENDDO 


DO IZ = l Mz 


CSPD (12,1, IZ) = AVES (IZ) 
CSPD (11,2, IZ) "=A AVES (IZ) 
CSPD (10,3, IZ) = AVES ra) 
CSPD(9,4,IZ) = AVE5 (IZ) 
CSPD (8,5,IZ) = AVE5 (IZ) 
CSPD (7726, στ) ταν τὸ) 
CSPD(6,7,IZ2) » AVES(IZ) 
CSPD(5,98,IZ2) = AVES( 12} 


CSPD(4,9,I2) = AVES (12) 
ՇՏՔ. 1210, 12) = AVES GZ) 


CSPDO2,Il1512) -JAVESTPZ) 
CSPD( ”1, 12; 12)mm AVES (ez) 
CSPD(13,1,IZ) = AVE4(IZ) 
CSPD(12,2,IZ) = ՃՄՔ4(12) 
С5РО(11,3,12) = AVEB4 στ. 
CSPD(10,4,IZ2) = AVE4 (IZ) 
CSPD(9, 5,12) = AVE4( 22) 
CSPD(876,12Z)) = AVE4 (i> 
CSPD (7) 1. 2) = ΑΕ. το 
CSPD(6,8, LZ) = AVE4 (12) 
CSPD(5, 9. 12) = AVE 3 (2) 


CSPD(4, 10,12) ΑΕ. πο 
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CSPD(3,11,IZ2) = AVE4 (IZ) 
ESDDI(2,12, IZ) = AVE4 (IZ) 
ESDD(I,13,IZ) € AVEA(IZ) 
CSPD (13,2, IZ) = AVE3 (IZ) 
ЕСВ0012,3,12) = АУЕЗ (І2) 
CSPD(11,4,IZ) = АУЕЗ (12) 
CSPD (10,5, IZ) = AVE3 (IZ) 
CSPD(9,6,IZ) = AVE3 (IZ) 
CSPD(8,7,IZ) = AVE3(IZ) 
CSPD(7,8, IZ) = AVE3(IZ) 
CSPD (6,9, IZ) = AVE3(IZ) 
Sopp (5,10,1Z) = AVE3(IZ) 
ՇՏՔը (4,11,12) = AVE3(IZ) 
ΟΡ 12,17) = AVE3 (IZ) 
CSPD(2,13,IZ2) - AVE3(IZ) 
ESPD(I3, 3, IZ) = C2 (IZ) 
BSDD(I2,4,IZ2) = Շ2(12) 
ESPD(ILI,5,IZ) * C2(IZ2) 
CSPD(10,6,I2) - C2(IZ) 

S PD(9 7,17) = C2(I2) 
ESDD(8,8,I2) - C2(IZ) 
ESDO(C,9,I2) *-cC2(IZ) 
CSPD(6,10,I2) = Շ2(12) 
CSPD(5,11,12) - C2(IZ) 
CSPD(4,12, IZ) = C2(I1I2) 
EC PD(3,13, IZ) = C2(I2) 
SPD(13,4, IZ) = C2 (I2) 
ESPD(12,5,I2) - C2(I2) 
ESPD(II,6,IZ) 3 C2(IZ) 
ОРО (10,7,ї2) - C2(IZ) 
CSPD (9,8, IZ) = C2(IZ) 
CSPD(8,9,IZ2) - C2(IZ) 
С5РО(7,10,2) = ԸՇ2 (12) 
C PD(6,11,172) =.C2 (IZ) 
С5РО(5,12,12) = С2(17) 
ШОЕР(1,13,т2) = C2(IZ) 
ШОРО(13,5,т2) „= С2(т2) 
ESED(I2,6,IZ2) = C2(IZ) 
ШОРЫ 11,7,12) = С2(т2) 
CSPD(10,8,12) - C2(IZ) 
ՇՏՔը (9,9,12) » Շ2(12) 
ՇՏՔը (8,10,12) - C2(I2) 

— CSPD(7,11,IZ) - C2(IZ) 
ՇՏՔը (6,12,12) Հ C2(IZ) 
τ! 13,172) = C2(IZ) 
CSPD(13,6,1Z) = AVE2(IZ) 
СЅРр (12,7,12) = АУЕ2 (12) 
ШБРр(11,8,Т2) = AVE2 (IZ) 
CSPD(10,9,IZ2) - AVE2(IZ) 
CSPD (9,10, EZ) = AVE2 (IZ) 
CSPD (8,11, IZ) = AVE2 (IZ) 
CSPD (7,12, IZ) = AVE2 (IZ) 
CSPD (6,13, IZ) = AVE2 (IZ) 
CSPD(13,7,1Z) = AVE1 (IZ) 
CSPD(12,8,1Z) = AVE1 (IZ) 
CSPD(11,9, IZ) = AVE1(I2) 
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CSPD(10, 10. D2) = AVE ЕУ) 


CSPD(9,I1,IZ) = АЕ (ТА) 
ՇՏՔՇ (8, 12, 12) » ՃՄՔ1(12) 
ՇՏՔը (7,13, 12) Հ ՃՄՔ1(12) 
CSPD (13 39 712) == C τα, 
CGCSPD( (12,912) = сүт?) 
ο ο 10, ΤΘ] 61 (17) 


ΟΣΕΡ 11, 125 851112) 
CSPD09: 12212) "ΕΠΙ ΘΙ 
CSBDI(S, 13, IZ) = C1 (12) 
CSPD (13,9, I2) =m. C1 (I2) 

CSPD(12/10,12)>”= CEA) 
CSPD (11, 11, T2) Se (T2) 
CSPD (10,12, 12) C LT) 
CSPD( (9,13, 12) = 61112) 


С5РЮ (13,10,12) ο τα 
ՇՏ (12, 11, Լ ՀԱԱ ՏՆԵՆ (ТА ) 
CSPD (11; 12 ΤΩ) ΞΕΡΕΙ 
CSPD (10) 137.12) s h ro) 
CSPD (13 111225 τα 
CSPD(12, 12,22) ο τα) 
CSPD (11,13,12) + C1(IZ) 
СРБ (1З, Т2, ТЈ АИСТ ТА) 
CSPD/(12, 13. 12) ο τα) 
CSPD(13,13, 12) &= 21 (12) 


ENDDO 


k k k k k k k k k k k kk k k k k k k k k k k í k k k k k k k k k £ k k £ k £ k k k í £ £ k k k k k k k k k k k k £ k k k k k £ £ k k k k k k k k k 


C WRITTING TO THE INPUT FILE FOR GRIDDER. (GRIDSVP.DAT) 


k k k eoe oe oe eoe ve € ΧΧΧ 6 ke oe eoe oe € ok € e ke oko e ee oe oe e eoe eoe € e ok n xxx xx kx x Χ 


WRITE(74) (XCOORD(I), I-21,MX) 
WRITE (74) (YCOORD (J), J=1,MY) 
WRITE (74) (2С(К), K=1,MZ) 
Ρο =. 1, Mx 

DO J = 1,MY 

WEITE(T4) (CSEDUDL Τ.Ε), Κι м2) 
- write(6; 7) թը (T J h) 

ENDDO 

ENDDO 


X koe o e ecco ee eoe Kx xx XK xxx xx xx xx xx xxx xxx xx xx xx |J x € kx kx kx kx AG x kx kx xxx xxx xxx xxx xxx 


CLOSE (41) 
CLOSE (42) 
CLOSE (43) 
CLOSE (44) 
CLOSE (74) 


STOR 
END 
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hyd8beam.m 


ο9 


Hydrophone #8 Arrival Plot 


* k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k f k k k kf k k kk k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k Kk 


ա Ա: J.M. Elliott, LCDR/USN, NPS, 1992 


Purpose: 

This Matlab program plots HARPO ray arrival structure at a particular 
receiver hydrophone array location. Plotting time vs relative amplitude for 
the eigenrays arriving within one wavelength of the hyd depth and from 
a specified arrival angle. Rays that follow essentially the same path are 
taken as one (this is indicated by rays with the same # of turning points, 
surface bounces, and bottom bounces). 

Two different plots are produced: 

1) representing received pulses as unit impulse and 
2) representing receiced pulses as gaussian distributed in time with 
a pulse width equal to 1/(band width). 

Input is data in f3res.dat which is output from the Fortran program 

cordat.f which is corrected HARPO vertical recvr surfaces data. 


F— o9 o9 o9 o9 o9 o9 oO o9 o9 oO oO g9 g6 o9 οἱ οἱ οἱ оо 


oad f3res.m $ f3Jres.m -» f3res.dat 
Դ րրթ- ths, thr, z, tpcnt, srfcnt, botcnt, ramp, t 


Find the eigenrays arriving within the desired window around the hydrophone. 
Compute the relative amplitude and total TL computed by cordat.f and plot the 
data: 

Define: 

Z = arrival depth of eigenrays in km above mean sea level (msl) (а neg ἢ) 

b = index vector corresponding to all rays arriving w/i 1/2 wavelength 
of the hydrophone 

c = index vector corresponding to the rays arriving w/i 1/2 wavelength 
of the hydrophone and coming in from above (neg grazing angle - 
<-0.00009) or below (pos grazing angle - >0.00009) the horizontal. 

t = eigenray arrival time in sec 

D = pulse width (sec) = 1/band width (hz) 

ο. 1/16, 

tlt = total trans loss in dB (absorption is neglected) 


ОО οἱ оо o9 o9 оо o9 o9 o9 g9 o9 оо 


oo οὐ 





sls = signal level lm from the source in dB 
sls - 183; 

5 slr = signal level at receiver in dB 

% amp - amplitude of eigenrays at receiver in Pascal 

% s = amplitude of received signal simulated as gaussian pulses 

Ф 2 = Һуагорһопе depth above msl in km (a negative number) 
0.220; 

% lambda = acoustic wavelength in km = (c(average)/f)*1e-03 
lambda = (1445/224) *le-03; 

% uz = upper limit of window 

| uz = z + lambda/2; 
5 lz = lower limit of window 


12 = z - lambda/?; 


Seer ores (:, 3); 

find(z<=uz & z>=l1z); 

@eomeeing (f£3res (b(:,1),2)>0.00009) ; 
n,m} = size(b); 


e QU O N 
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Define a new matrix "resd(n,7)" for ease in displaying the output data 
for hydrophone. 
resd = ones(n,7); 


5 
5 


resd(:,1) = f3res(b(:,1),1)*resd(1,1); $ launch angle (deg) 
resd(:,2) = f3res(b(:,1),2)*xresd(1,2); % arrival angle (deg) 
resd(:,3) = f3res(b(:,1),4)*resd(1,3); Ф по of turning pomikss 
resd(:,4) = f3res(b(:,1),5)*resd(1,4); Ἔ no of surface bounces 
resd(:,5) = f3res(b(:,1),6)*resd(1,5); % no of bottom bounces 
resd(:,6) = f3res(b(:,1),7)*resd(1, 6); $ relative amplitude (dB) 
resd(:,7) = f3res(b(:,1),8)*resd(1,7); $ arrival time (sec) 


% Now sort the rays in order of increasing arrival time: 
[Y,I] = sort(resd(:,7)); 


EIGENRAY DETERMINATION: 

Eliminate duplicate launch rays and rays that are essentially the same 
due to identical number of turning points, surface bounces and bottom 
bounces. First define a new matrix resn with the data in order of 
arrival time: 

resn = ones(n,/7); 

resn(:,1)=resd(I,1) *resn(1,1); 
resn(:,2)=resd(I,2) *resn(1, 2); 
resn(:,3)=resd(I,3) *resn(1, 3); 
resn(:,4)=resd(I,4) *resn(1, 4); 
resn(:,5)=resd(I,5)*resn(1,5); 
resn(:,6)=resd(I,6) *resn(1, 6); 
геѕп(:,7) =геѕа(ї, 7) *геѕп(1, 7); 


o? o9 


a of оо 


for i τα; 
for k = 1:n; 
1Ё 1 == k; 
else; 
difa = abs(resn(i,l) - resn(k,1)); 
diftp = abs(resn(i,3) - resn(k,3)); 
difs = abs(resn(i,4) - resn(k,4)); 


difb = abs(resn(i,5) - resn(k,5)); 
1f ((difa == 0.0000) | ((diftp==0) & (difs==0) & (difb==0))); 
; resn(i,:)=({0 00000 Oj; break; end; 


end; 
end; 
end; 
Now eliminate the unwanted ray zero data lines and call new data 
matrix resf: 
К = 1; 
for i = ո: 
if resn(i,1)~=0; 
resf(k,:)=resn(i,:); 
к = BUS 
епа; 
епа; 


oo 


oe 


5 Now plot the eigenray arrival time vs relative amplitude: 
[n,m] = size(resf); 
C= οποια N 


πο 
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tlt = ones(n,1l); 
slr = опеѕ (п,1); 
amp = опеѕ (п, 1); 
mor: = 1:n; 
muy) = (resf(i, 7) *t(i,1)); 
% mime!) = (f3res(b(c(i,1),1),5))*tlt(i,1):; 
% СИЕЗ l) - (sls - tlt(i,1))*slr(i,1); 
amp(i,1) = ((1le-06) *10% (resf(i,6)/20))*amp(i,1); 
$ if slr (i,1l1)>=100; 
% hea, |) 07 
$ end; 
end; 


бахіѕ (134.62 36.0 0 0. 01]); 
ESEEpISGt(211),plot(t,amp,'.*') 
title('HYD #8 EIGENRAY ARRIVALS') 
ο τς 4.65, «011, 517’) 

ο 5.090.007, 3185”) 
8:535. 32,0. 0085, 7) 
' ξεκε(35.80,0.003,’Β20’) 
| xlabel('Arrival Time (sec)') 

| ylabel('Relative Amplitude (Pa)’) 
| $text (36.5,2,'Hydrophone Depth = 220m’) 
$text(36.5,-2,'AÀArriving w/i 1/2 Lambda') 
| %šaxis([1 2 3 4]);axis; 
$pause; 
51212), բՔ1Շե (Էէ, Օ6Տք(:,6),՛.") 
meext (34.65,82,’R17’) 
meext (35.03,78,'R18°) 
8500055. 32,80," R19') 
8e xE(35.80, 70, ' R20^) 
xlabel('Arrival Time (sec)’) 
|ylabel('Relative Amplitude (dB)’) 


Simulate received eigenrays as gaussian pulses: 
Define: 
ta = time axis 
gisa = gaussian intermediate signal amplitude 
gsa = gaussian signal amplitude 
mas. O00: .002:37.5000; 
(o,p]=size (ta); 
gisa = ones(n,p); 
gsa = ones(1,p); 
Е 
а 


o9 o9 оо оо 


c? o9 o? o? 


ΙΙ 7.40 35.0917 35.3814 35.8676]’; 
fee | Le-0'6*10°(79.5235/20)  165-06»10՞ (75.6184/20).. 





% meso lO 77. 3569/20) 16-06710 (66.5519/20) ] 
.- ) = 1:р; 
ա 5 = 1:n; 
% if і== 
Е EE cnn) -exp(—2~((ta(l, })—-t (1,1))*2)/ (D*2))*gisa(i, 4): 
7 end; 
{ 1-1: 
E Ա ԱՆ-Ն d *expi-2*titat1,3)-E(1,1))^2)/(D^2)) *gisa(i,j)... 
| 5 αρα ιτ ο) 
1$ end; 


[5 end; 


po 
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$end; 

$gsa = sum(gisa); 

$axis([34 255 ο ο αἱ]. 
$subplot(212) 9 )pilot(ta,gsa) 
$title('SIMULATED PULSE ARRIVALS’) 
$xlabel('Arrival Time (sec)') 
$ylabel('Relative Amplitde (Pa)’) 
$text (36.3,1e+04,’Pulse Width = 62.5 msec’) 
$axis([1 2 3 4));axis; 

Stext (36,-68,’Launch Angles 0-25 Deg’) 
Stext (36,-71,’Launch Step = 0.01 deg’) 
$text (36,-74,’Broad Front (40km)’) 
$text (36,-77,’Hyd Depth = 150m’ ) 


74 


SN - (fft(cn,512))/sum(ac); 
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Qe 


beampat.m 


% Barents Array Beam Pattern 

Qe k * k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k kk k k k k k k k k k k k k k k kk kk k k Kk k k kk xxx k 
% Written by: J.M. Elliott, LCDR/USN, NPS, 1992 

% Purpose : 

5 This Matlab program plots an array plane wave beam pattern based on 

$ a method for a linear array of point hydrophones by Ziomek (1985). 

% 

go * K k k k k k k k ` k k k k k k k e oe ЖК 
N = 16; $ М = # of hydrophones 

d = 10; % d = hydrophones interelement spacing (m) 
ք = 224; Ф Ё = frequency of interest (Hz) 

с = 1450; % с = sound velocity (m/s) 

for steer = 20 ; % steer = angle in deg main lobe is to be 


steered off broadside of array 
Sip = steered angle ref to broadside 
xp = steered angle in ‘u’ space 


psip = 90 + steer; 

fxp = f*cos (psip*pi/180) /c; 
ñ = —(N-1)/2:1: (N-1)/2; hyd # index 

ac = cos ( (cos ( (pà* ( (2*n-1)/2) /N)))): $ ac = cosine amplitude window 
theta = (-2*pi*fxp* ((2*n-1) /2)*d); % theta = phase weight 
% 
% 


оо оо ое 
J r'U 


Сп =.ac’ .*exp(j*theta’); cn = complex wieghting 

SN = normalized directivity 
% function 

БИО 1:511; t$ me«sFFT bin +t 


m 
ΠΠ. (acos((c.*m')/(£*512*d4)))*(180/p1); $ psi = conversion fm 'u' to deg 
$ nw otc.*m')/(f*224*d); 


end; 

plot(psi,abs(SN)) 

title('PLANE WAVE BEAM PATTERN') 
xlabel(’Degrees From Axial (90=broadside, QO=endfire)’) 
ylabel('Normalized Amplitude’ ) 

text (75,0.9,'16 Element Vertical Linear Array’) 
text (75,0.85,’10m Interelement Spacing’ ) 
Eext(75,0.8,' £ = 400 ΗΖ’) | 
text(75,0.75,'Cosine Amplitude Window’ ) 
$text(73,0.5,'Steered to 20 deg off Broadside’ ) 
text(75,0.7,' Beam Width = 1 deg’) 

text (80,0.3,’Grating Lobe 21 deg away’) 

ѕепа 


15 
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resanl.m 


Resovability Analysis Table Generator 


* e oe k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k x xA k k k | 


Written by: J.M. Elliott, LCDR/USN, NPS, 1992 


Purpose: 

This Matlab program uses the output data file <f3res.dat> from the 
fortran program <cordat.f> and provides an output data table of tomograph- 
ically resolvable rays for the whole array as a plane wave beamformer. 

To compute the data table of tomographically resolvable rays for each 
hydrophone as an omnidirectional receiver, the program should be modified 
as indicated in the comment lines in the body of the program. 

The data is presented in the fol format columns: launch angle, arrival 
angle, no of turning pts, no srf, bounces, no bot bounces, rel amp and 
arrival time. 

The data is displayed in order of increasing arrival time and only for 
rays arriving w/i 1/2 wavelength (of the freq of interest) from each 
hydrophone. 

Rays following essentially the same paths, as indicated by similar 
number of turning points, surface bounces, and bottom bounces, are taken as 
опе ray (the first ray of the group arriving in time). 

Each ray is tomographically resolved in time using a specified pulse 
width. 


wae kk kee kk kh kkk kk kk kkk kkk kkk kkk kkk kkk kkk kkk kkk kkk Kk Kk kh kkk KKK kkk KKK KKK KKK KKK K 


Define: 
z = hydrophone depth above msl in km (a negative number) 
lambda = 1/2 acoustic wavelength in km = 0.5 * (c(average)/f)*1e-03 
lambda = 0.5*(1445/224) *le-03; 
uz = upper limit of depth window 
uz = z + lambda; 
lz = lower limit of depth window 


% 

q 

% lz = 2 - lambda; 

% ds = Arrival angle (grazing) sector. Used for plane wave beam- 
$ formed array table only. One of 15 defined as follows: 
$ (Three degree beam widths used based on 224 Hz plane wave beam 
5 pattern analysis) 

Հ ՅՑ 122.56 ----” 

% = d2 (-19.5 to -16.5) 

% = аз (πας ο б =13 5) 

% κ οι (213.5 tor ME 

$ -cd5- (095 το. ο) 

% = d6 (-7.5 to -4.5) 

% = d7 (-4.5 to -1.5) 

5 ια, кк ә е е с 

% = d9 (+1.5 to +4.5) 

% = ТО (44.5 τσ ο. 

> = dll το το πο) 

2 ա զ12 (10 50060 Են 5 

7 ՏՎ13 το δι 755 

? = 014 πι ο τσι ο 5) 

7 -α15.{13 5 Έα. 505) 

load f3res.m $ f3res.m - f3res.dat 
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ջ,»թՉ29--Լա ths, thr, z, tpcnt, srfcnt, botcnt, ramp, t 

$ Let z = the third column arrival depths: 

,-ծ.`. 3-65 (:, 3): 

$ 

ls Determine all eigenrays or all rays arriving w/i 1/2 wavelength of the 
|$ frequency of interest in depth of the 16 hydrophones: 

М First define hydrophone depths in km above mean sea level: 

oi = -.150; 





122 = -.160; 
ԱՅ - -.170: 
z4 = -.180; 
25 = -.190; 
26 = -.200; 
ЕГ = -,210: 
mo = -.220; 
29 = -.230; 
‚210 = -.240; 
mil - —.250; 
212 = -.260; 
Ё13 = —.270; 
zl4 - -.280; 
215 = -.290; 
216 = -.300; 


$ Next define lambda depth window around each hydrophone: 


uzl - zl Հ lambda; 
Ш21 = 71 — lambda; 
022 = z2 + lambda; 
1z2 = z2 - lambda; 
uz3 = 23 + lambda; 
123 = z3 - lambda; 
024 = z4 + lambda; 
1z4 = 24 - lambda; 
uz5 = z5 + lambda; 
125 = 25 - lambda; 
026 = 26 + lambda; 
126 = 26 - lambda; 
027 = z7 + lambda; 
127 = 27 - lambda; 
028 = 28 + lambda; 
128 - z8 - lambda; 
ս29 - 29 + lambda; 
129 - z9 - lambda; 
ս210 - 210 + lambda; 
1z10 = 210 - lambda; 
ս211 = z11 + lambda; 
1211 = 211 - lambda: 
B212 2 212 + lambda: 
1212 - z12 - lambda: 
0213 = zl13 + lambda; 
1z13 = 213 - lambda; 
ս214 = 214 + lambda; 
1z14 = z14 - lambda; 
ս215 = 215 + lambda; 
1215 = 215 - lambda; 
uz1l6 = 216 + lambda; 


πο 
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1216 = 216 - lambda; 


% Now find the matrix index number associated with the eigenrays arriving 
% w/i the depth windows around each hydrophone. For individual hydrophone 
$ table generation all but the hydrophone of interest must be commented out 
$ and the program run individually for each hyd, exiting and reentering 

$ Matlab each time. 

bl*iEBgInd(z«uzi & z>121 


b2 = find(2<022 &.Z2?122), 
b3 = find(2<u23 & 22123); 
b4 = find(z<uüu2z24 6 22124): 
δα Επ νο οἱ 
Ьб - find(z«uz6 & z»l1z6): 
b7 = find(z<uz?7 & 2>177); 
b8 = find(z<uz8 & z>1z8); 
b9 = find(z<uz9 & z>1z9); 
b10 = find(z<üz10 & 2>1z10); 
b11 = find(z<üuz1ll & 251211) 
b12 = find(2<u212 4 2> IE 
b13 = епа оз КИРЕ 
Ե14 2 ք1ոՎ(2ՀՀս214 & 251214): 
ԵԼ5աքոՎմ(22ս215852-1215)7 
b16 = findad(z<üz16 0802 1215); 


Ф Let matrix c = all rays arriving w/i 1/2 wavelength of all 16 hyd’s. Comment 
$ out this step when computing individual hyd resolvability table. 
= (bl; b2; b3; b4; b5; b6;-b?; 58; ΕΕ ΕΙ, Ա IE 

ББТ 


Q 


Now determine the rays arriving in which of the 15 beam sectors of 
three degree beam widths each. The program must be run individually 
for each sector and Matlab must be exited and reentered each time. 
Comment out this step when computing individual hyd resolvability table. 
- find(f3res(c(:,1),2) »2-22.,5000.8$ £3res(c(:,1)92) «19 ο. 
find(f3Jres(c(:,1),2)»-2-19.5000 & f3resí(c(:,1),2)«-16.5000DEE 
find(f3res(c(:,1),2)»-—-16.5000 & f3res(c(:,1),2)<=13 S000- 
find(f3res(c(:,1),2)>=-13.5000 & f3res(c(:,1),2)<-10.5000); 
find(f3res(c(:,1),2)»-2-10.5000 & f3res(c l; ο s sO 
find(f3res(c(:,1),2)>=-7.5000 & f3res (c(:/1),2)<°430n5000) 
find(f3res(c(:,1),2)>=-4.5000 & f3res(c(:,1),2)<-1.5000); 
find(f3res(c(:,1),2)>=-1.5000 & f3res(c(:,1)y2)<t105000) 
find(f3res(c(:,1),2)>=+1.5000 & f3res(c(:,1),2)<+4.5000); 

find(f3res(c(:,1),2)>=+4.5000 & f3res(c(:,1)” 2) οσα, 
find(f3res(c(:,1),2)>=+7.5000 & f3res(c(:,1),2)<t10 5060005 
find(f3res(c(:,1),2)>=+10.5000 & f3res(c(:,1);2)S 115 ΠΠ. 
find(f3res(c(:,1),2)>=+13.5000 & f3res(c(:,1),2)<+16.5000); 
find(f3res(c(:,1),2)>=+16.5000 & f3res(c(:,1),2)<+19 5070D 
find(f3res(c(:,1),2)>=t19.5000 & f3res({(c(: ,1)72)<-22 50. 


l #H H Hl H H H MW 


oO oO oO օօ 30 30 оо оо оо SQ oQ oO οἱ SQ (L o9 o9 o9 oO 


| 

az 
аз 
d4 
d5 
d6 
а7 
а8 
а9 
dı 
dl 
di 
dl 
dl 
al 


0 
1 
2 
3 
4 
5 


Define a new matrix "resd(n,7)" for ease in displaying the output data 
for sector dS. Sector dS must be changed for each i programs 
When computing individual hyd resolvability sector dS should be 
set equal to the hydrophone of interest, bX, vice sector, dx. 
аз = ат 
[пуп] = втте(а1):; 


oP of alO 39 
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геѕа = опеѕ (п, 7); 

--շաՆ 1 թյ) տ f3res(c(dS(:,1),1),1)*resd(1,1):; $ launch angle (deg) 
$ for hydrophone resolvability table generation use the fol format: 
,Խ,--Հ(:, 4)" £3res(dS(:,1),Y) *resd(1,X): 


E0772)» f3res(c(dS(:,1),1),2) *resd(1,2):; $ arrival angle (deg) 
resd(:,3) = f3res(c(dS(:,1),1),4) *resd(1,3): % no of turning points 
resd(:,4) » £3res(c(dS(:,1),1),5) *resd(1,4); $ no of surface bounces 
nesd( 5) 9 Ff3res(c(dS(:,1),1),6) *resd(1,5):; $ no of bottom bounces 
resd(:,6) = f3res(c(dS(:,1),1),7)*resd(1,6); Ф relative amplitude (dB) 
resd(:,7) = f3res(c(dS(:,1),1),8) *resd(1,7); $ arrival time (sec) 


Ф Now sort the rays in order of increasing arrival time: 
[Ү,Т] = sort (resd(:,7)); 


$ Eliminate duplicate launch rays and rays that are essentially the same 
| $ due to identical number of turning points, surface bounces and bottom 
% bounces. First define a new matrix resn with the data in order of 
$ arrival time: 

resn - ones (n,7); 

resn(:,1)resd(I,1)*resn(1,1); 

resn(:,2)=resd(I,2) *resn(1,2); 

resn(:,3)-resd(I,3)*resn(1,3): 

resn(:,4)=resd(I,4) *resn(1, 4); 

resn(:,5)=#resd(I,5) *resn(1,5); 

resn(:,6)=resd(I, 6) *resn(1, 6); 

resn(:,7)-resd(I,7)*resn(1,7); 

ου. του]. 

Or k ż l:n; 


if i == k; 
else; 
difa = abs(resn(í,l) - resn(k,1)); 
diftp = abs(resn(i,3) - resn(k,3)): 
difs = abs(resn(i,4) - resn(k,4)); 


difb = abs(resn(i,5) - resn(k,5)); 
2£ ((difa <= 0.0000) | ((diftp<=0) & (difs<=0) & (difb<=0))); 
; resn(i,:)=({0 0 00 00 Oj; break; end; 


end; 
end; 
end; 
Now eliminate the unwanted ray data lines: 
k = 1; 
hor 1  l:n; 
if resn(i,1)-20; 
resf (k,:)=resn(i,:); 
k = k + 1; 
end; 
end; 


оо 


Now determine tomographic temporal resolvability by comparing each 
ray arrival time to its neighbor for a distance of the pulse width of 
interest. If the distance is »- the pulse width then the ray is 
resolvable (R), if not the ray is not resolvable (U). 

224 Hz source - Pulse width 1/ (band width) 
400 Hz source = Pulse width 1/ (band width) 


o9 o? 


o9 o9 


1/16Hz = 0.0625s 
1/100Rz = 0.010s 


o9 o9 
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The new [n,1] matrix resg contains the 'U' for unresolved or 'R' for resolved 
info corresponding to the n rows of the resf matrix. 
[n,m] = size(resf); 
fons: = lean: 
for k= Wen. 
if i ==k; 
else; 
dift = abs( resf(i,7)-resf(k,7)); 
if dift <= 0.01; resg(1,1)='U'!'; break; end; 
if dift > 0.01; resg(i,1)='R'; end; 
end; 
end; 
end; 


ae oP 


To get a print out of the plane wave beamformed array resolvability data 
table for the sector of interest (or of the individual hydrophone data 
table for the hydrophone of interest) sorted by time type: resf. 

To get a print out of corresponding row 'U' or 'R' info type resg. 


oO oO oO oP 


80 


APPENDIX B: INDIVIDUAL HYDROPHONE RESOLVABILITY ANALYSIS 
TABLE 


The following table was generated by treating each 
hydrophone of the 16 element vertical array (B) as an 
independent omnidirectional receiver. The table was compiled 
using the MATLAB computer program касаба: А copy of 
'resanl.m' can be found in Appendix A. Definitions of RAMP, 
R, and hydrophone depths can be found at the end of the 


Appendix. 
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LAUNCH 
ANGLE 
(deg) 


ա ա C29 Cà CO CO CO CO Www Ow O1~) © AAO OV 


КӘ tS ԷՉ Кә թ տ տա (л (л O © 


„9401 
„8401 
‚1600 
„1500 
„3300 
„3601 
.2100 
.7400 
.0600 
.5200 
‚5100 
„0300 
‚0700 
‚0200 
„4100 
. 7300 
.1100 
„1200 
‚1400 
‚4900 
. 1401 
.8301 
.4301 
.4500 
.5900 
.6001 
„4401 
„8700 
‚1901 
„9101 


‚8801 
.4201 
.2600 
.4100 
.7800 
.0800 
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HYDROPHONE 11 


10.2900 28955201 22.0000 4.0000 13.0000 27495918 34.5544 
7.2901 6-9111 19.0000 0 10.0000 176535183 34.5597 
NOH Z95229 24.0000 4.0000 13.0000 74.1964 34.5874 
4.1100 4.2339 32.0000 0 13.0000 1517570 34.6069 
0.1200 2.2687 28.0000 0 17.0000 87.5601 34.6085 
9.4101 8. 6182 15.0000 8.0000 12.0000 74.1037 34.6112 
2.6100 =3.1218 33.0000 0 14.0000 81.5199 34.6118 
2.2900 2.8892 36.0000 0 16.0000 77.2450 34.6125 
1.5100 2.6044 32.0000 0 19.0000 84.6509 34.6129 
1.9700 2.8663 31.0000 0 17.0000 88.7595 34.6145 
1.9800 2 58316 32.0000 0 17.0000 86.3594 34.6145 
9.2601 8.6170 16.0000 8.0000 12.0000 71.0924 34.6289 

12.5001 11.1888 0 21.0000 21.0000 65.6602 34.9902 

11.8600 -=11 3191 0 21.0000 20.0000 71.5419 35.0498 

14.0601 14.0723 0 26.0000 26.0000 78.8534 35.4251 R 

13.1501 114151959 0 26.0000 25.0000 121123 35.5364 К 

16.7200 16.5431 0 32.0000 32.0000 66.8359 35.9145 к 

15.9200 -16.5941 0 32.0000 31.0000 63.4517 36.0304 R 

24.6901 20.3055 0 44.0000 44.0000 14.0856 36.3099 

23.3001 20.1363 0 42.0000 42.0000 35.7012 36.3446 

22.7001 -20.2556 0 42.0000 41.0000 35.6383 36.4617 R 

20.3200 19.7987 0 41.0000 41.0000 39.5796 36.8971 R 

19.6601 =-19. 9250 0 41.0000 40.0000 40.5354 37 012885 

HYDROPHONE 12 
6.4301 6.2827 24.0000 0 9.0000 89.6127 3455] 
5.8800 5.8484 41.0000 0 12.0000 79.4809 34.6035 
4.4000 4.3367 32.0000 0 13.0000 81.1690 34.6044 
5.1800 -5.0945 34.0000 0 12.0000 77.9052 34.6076. 
0.1500 2:2201 27.0000 0 17.0000 88.5312 34.6085 
0.1600 2.1676 . 28.0000 0 17.0000 86.2914 34.6085 
0.1700 Հան թել 29.0000 0 17.0000 87.1690 34.6085 
3.3200 23:6215 34.0000 0 13.0000 80.8535 34.6087 
3.2400 3.5394 32.0000 0 14.0000 79.4335 34.6113 
2.6000 =2.9887 3620000 0 14.0000 79.7294 34.6118 
1.4800 295083 31.0000 0 18.0000 86.1824 34.6128 
2.1300 -2.7349 34.0000 0 16.0000 81.7180 34.6134 
2.0000 2081537 34.0000 0 17.0000 83.0568 34.6147 
7.7101 =8 3325 26.0000 0 12.0000 79.3245 34.6980 R 

11.9700.  —P52825 0 21.0000 20.0000 7879192 35.0374 

11.8300 11. 1377 0 21.0000 21.0000 70.3082 35.0798 

13.9401 14.0684 0 26.0000 26.0000 78.8688 35.4415 R 

13.1801 =14 1933 0 26.0000 25.0000 1255135 35.5344 R 

16.6500 1655 52 0 32.0000 32.0000 66.1088 3579256845 

15.9700 -1675976 0 32.0000 31.0000 659-9213 ος ο κ 

24.5101 20.2899 0 43.0000 43.0000 26.8543 36.1844 R 

23.2701 20.1347 0 42.0000 42.0000 35.6155 36.3501 R 

22.1301 —20.2143 0 42.0000 41.0000 35. 7554 36.4572 R 

20.2800 19.7971 0 41.0000 41.0000 39.4713 36.9041 R 

88 
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19.7501 -19.8870 0 41.0000 40.0000 41.1746 36.9950 R 
HYDROPHONE 13 
7.3201 6.8537 20.0000 10.0000 75.4929 34.5577 
9.8301 8.5903 16.0000 9.0000 13.0000 72.4467 34.6068 
4.1200 4.1334 33.0000 0 13.0000 79.4189 34.6073 
4.1300 4.1085 31.0000 O 13.0000 84.1847 34.6076 
| 0.2100 2.0172 27.0000 0 17.0000 87.7163 34.6083 
| 0.9000 -2.1376 32.0000 0 16.0000 92.4143 34.6090 
| 0.9800 -2.1914 30.0000 O 16.0000 93.9919 34.6091 
| 1.0100 -2.2061 28.0000 О 17.0000 96.1761 34.6093 
| 1.0200  -2.2120 29.0000 0 17.0000 93.2634 34.6093 
| 1.0400 -2.2169 31.0000 0 16.0000 95.1576 34.6093 
| 1.0500 -2.2217 30.0000 O 17.0000 92.2944 34.6093 
| 3.3000 -3.5531 32.0000 О 13.0000 83.7259 34.6093 
| 1.0700 -2.1951 27.0000 О 16.0000 90.1986 34.6094 
| 1.0800 -2.1601 29.0000 O 16.0000 88.1058 34.6095 
| 2.3000 - 2.6582 34.0000 0 16.0000 77.3741 34.6121 
|. 1.4500 2.3618 31.0000 0 18.0000 78.9123 34.6125 
| 2.1200 -2.6580 36.0000 0 16.0000 81.8799 34.6132 
2.0100 2.6206 34.0000 О 17.0000 80.8960 34.6145 
9.2101 -8.4851 22.0000 3.0000 12.0000 85.5376 34.6262 
7.7201 -8.3059 26.0000 12.0000 83.3586 34.6982 R 
12.0100 -11.2540 0 21.0000 20.0000 70.3709 35.0335 
12.0300 11.1396 - Q 21.0000 21.0000 73.0061 35.0511 
13.8301 14.0650 0 26.0000 26.0000 78.6391 35.4564 Ք 
13.1901 -14.1535 0 26.0000 25.0000 70.6539 35.5337 R 
16.5700 16.5326 0 32.0000 32.0000 66.5419 35.9380 R 
16.0200 -16.5719 0 32.0000 31.0000 64.9706 36.0168 R 
24.5001 20.2913 O 43.0000 43.0000 26.1304 36.1849 R 
23.2401 20.1331 О 42.0000 42.0000 35.8631 36.3555 R 
22.7601 -20.2160 0 42.0000 41.0000 35.5000 36.4525 R 
20.2400 19.7955 О 41.0000 41.0000 39.8323 36.9108 R 
19.8001 -19.8891 0 41.0000 40.0000 41.0684 36.9859 R 
HYDROPHONE 14 
6.4801 6.2757 22.0000 0 9.0000 91.8753 34.5508 
6.2200 -6.1502 31.0000 0 8.0000 77.4682 34.5516 
6.4701 6.2828 46.0000 0 18.0000 81.7271 34.5520 
7.5701 -7.1387 22.0000 О 10.0000 70.7041 34.5776 
9.5101 8.5461 16.0000 8.0000 12.0000 70.1657 34.5985 
4.3800 4.2230 29.0000 О 13.0000 82.4682 34.6042 
4.1400 4.0386 32.0000 0 13.0000 82.6945 34.6073 
0.2600 1.9263 29.0000 О 17.0000 88.4186 34.6085 
0.8100 -1.9300 31.0000 O 16.0000 89.4463 34.6091 
0.8400 -2.0107 31.0000 O 17.0000 88.0979 34.6091 
1.1200 -2.0974 28.0000 0 16.0000 86.5873 34.6099 
1.1300 -2.0715 30.0000 O 16.0000 87.3547 34.6100 
1.4400 2.1630 30.0000 O 17.0000 83.6084 34.6120 
2.5300 2.8084 33.0000 0 15.0000 80.2181 34.6126 
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Оо 
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14.0000 
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7.0000 

8.0000 
12.0000 
12.0000 
17.0000 
17.0000 
16.0000 
16.0000 
17.0000 
13.0000 
13.0000 
16.0000 
17.0000 
17.0000 
13.0000 
16.0000 
12.0000 
12.0000 
20.0000 
26.0000 
25.0000 
32.0000 
31.0000 
42.0000 
41.0000 
41.0000 
40.0000 
39.0000 


16 
9.0000 


13.0000 
13.0000 


79 
82 
79 
70 
73 
85 
67 
78 
73 
66 
65 
22 
35 
36 
40 
40 


71 
72 
79 
78 
85 
88 
87 
88 
86 
71 
74 
83 
88 
83 
73 
80 
67 
80 
72 


13 
66 
64 
35 
J5 
40 
40 
33 


89 
69 
83 


.2481 
.0088 
.2159 
. 7126 
©3231 
. 7548 
. 3867 
.6462 
. 9508 
«4554 
«0836 
«8931 
«7261 
«71023 
«0015 
«8983 


«7002 
«2484 
«6993 
«6904 
«0701 
«9573 
«4888 
«0382 
«8456 
«71066 
«4812 
«4189 
«8305 
«9690 
«2620 
«6079 
«6847 
«93442 
«4460 
‚7124 
.4777 
. 3702 
.5200 
. 7648 
.8690 
.1053 
. 9451 
. 7800 


.8872 
. 2232 
.0594 


34 
34 
34 


.6130 
. 6132 
. 6141 
. 6421 
.0269 
.0636 
.1003 
.4699 
„5303 
. 9507 
.0097 
.1839 
. 3628 
.4444 
"9198 
. 9780 


.5146 
.5252 
.6001 
.6083 
.6084 
.6084 
.6090 
.6092 
. 6092 
.6093 
.6105 
«6107 
‚6118 
‚6118 
.6125 
.6129 
.6208 
.6987 
.0256 
.4819 
‚5202 
„9634 
«0033 
«4081 
«4395 
«9270 
«9690 
1354 


.9509 
«9540 
.6031 
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4.3100 3.9939 29.0000 0 13.0000 84.4584 34.6041 
4.2500 = 4. 0319 30.0000 0 12.0000 97.4540 34.6053 
4.1900 3.9263 31.0000 0 13.0000 84.1616 34.6054 
4.2400 at. 0008 28.0000 0 12.0000 83.1009 34.6064 
0.4000 155370 29.0000 0 17.0000 90.6496 34.6084 
0.4200 1.3062 30.0000 0 17.0000 9879022 34.6085 
0.7100 E7032 31.0000 0 16.0000 87.1150 34.6093 
0.7000 -1.6706 32.0000 0 16.0000 87.8697 34.6094 
3.2700 3.1659 36.0000 0 14.0000 81 .2573 34.6104 
1.2100 -1.9175 28.0000 0 16.0000 89.4330 34.6110 
NEU -1.9175 27.0000 0 16.0000 90.4993 34.6110 
1.2400 =1.8862 29.0000 0 16.0000 87.8468 34.6110 
1 3700 1.9304 33.0000 0 17.0000 86.8783 34.6119 
2.3200 2.3465 33.0000 0 16.0000 77.0698 34.6119 
J. 3500 1.8848 31.0000 0 17.0000 85 2 34.6120 
1.3600 1.9280 32.0000 0 17.0000 85.5294 34.6120 
255500 2.5488 30.0000 0 15.0000 78.9044 34.6122 
29 —2.3105 37.0000 0 16.0000 8075857 34.6131 
Ք -9.0491 17.0000 6.0000 10.0000 69 7367 34.6161 
9.2801 29:078 15.0000 8.0000 11.0000 5975166 34.6247 
արն -8.3734 28.0000 0 13.0000 74.1797 34.7021 R 
Ἱ 5501 14.0571 0 26.0000 26.0000 55 35.4936 
13.3401 -14.0889 0 26.0000 25.0000 74.2331 ου, οἱ 16 
16.3500 2525172 0 32.0000 32.0000 65.0820 35.9726 
' 0 -16.5484 0 32.0000 31.0000 64.9918 359978 
24.4701 -20.1854 0 43.0000 43.0000 23. 2387 36.1764 R 
s 570] -19.9860 0 42.0000 42.0000 35.0064 36.4325 R 
ΤΙ -19.6598 0 41.0000 41.0000 40.3247 36.9591 
րա Թ -19.8615 0 41.0000 40.0000 40.9984 36.9601 
ա ԷՈ -19.7090 0 40.0000 39.0000 36.1161 371353 R 


RAMP - relative amplitude in dB (RAMP=SL-TL, Sea State 3, 0.5 dB/bottom bounce) 


R - ray is temporally resolvable using a 62.5msec pulse separation. 


Hydrophone depths in km above mean sea level: 


,  Hydrophone f Depth 
al | 2.150; 

2 =. 160; 

3 -.170; 

4 -.180; 

5 -.190; 

6 -.200; 

7 -.210; 

8 -.220; 

9 2.2307 

10 -.240; 

11 -.250; 

| 12 -.260: 
T3 = 270: 

14 -.280; 

15 -.290; 

16 -.300; 


Note that the ocean bottom at the hydrophone array location is modeled at 318m 
Or -.318km above mean sea level. 


on 


APPENDIX C: LINE ARRAY RESOLVABILITY ANALYSIS TABLE 


The following table was generated by treating the 16 
hydrophones of the vertical linear array at position (B) in 
Figure 2 as a plane wave beamformer. The compilation was 
accomplished using the MATLAB computer program 'resanl.m'. A 
copy of 'resanl.m' can be found in Appendix A. Definitions of 
RAMP, R, and sector angles can be found at the end of the 


Appendix. 
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ARRIVAL 
LAUNCH ARRIVAL NUMBER NUMBER NUMBER TIME 
ANGLE ANGLE TURNING SURFACE BOTTOM RAMP (sec) 
(deg) (deg) POINTS BOUNCES BOUNCES (dB) (R-resovable) 
SECTOR 1 
EI 4701 -20.1854 0 43.0000 43.0000 23 2387 36.1764 
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MUI  -19.7114 0 40.0000 39.0000 33.7800 3771354 
SECTOR 2 
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S 9501 πο οσα ο 23.0000 3.0000 12.0000 2255930 ՅԾ ԵՈ 
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